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GLOBAL CLASSICAL SOLUTIONS TO PARTIALLY DISSIPATIVE 
HYPERBOLIC SYSTEMS VIOLATING THE KAWASHIMA CONDITION 

PENG QU AND YANJIN WANG 


Abstract. This paper considers the Cauchy problem for the quasilinear hyperbolic system of 
balance laws in R'*, d > 2. The system is partially dissipative in the sense that there is an 
eigen-family violating the Kawashima condition. By imposing certain supplementary degener¬ 
acy conditions with respect to the non-dissipative eigen-family, global unique smooth solutions 
near constant equilibria are constructed. The proof is based on the introduction of the par¬ 
tially normalized coordinates, a delicate structural analysis, a family of scaled energy estimates 
with minimum fractional derivative counts and a refined decay estimates of the dissipative 
components of the solution. 


1. Introduction 


Consider the following n-component quasilinear hyperbolic system of balance laws in d space 
dimensions: 


d 

dt{Giu))+J2{FHn))^ =S{u), 

k=l 

|f=0= Uq. 


( 1 . 1 ) 


Here u = (tti,..., : M x —>• R"' is the unknown, G = (Gi,...,Gn)^ : R” R”, 

= {F^, ..., F^f : R" ^ R" (A: = 1 , ..., d) and S = (5i,..., SnV : R" ^ R" are given 
smooth functions. It is assumed that VuG{u) is invertible, hence, by the chain rule, the system 
can be rewritten into 

d 

dtu + X ^^{u)uxk = Q{u), , . 

' k=i 

|t=0= Uo- 


Here the coefficient matrices A^{u) = iyuG{u))~^VuF^{u) {k = and the inhomo¬ 

geneous term Q{u) = {VuG{u))~^S{u). By hyperbolicity, for any u G R”' and any lo = 
(wi,..., UJd)'^ G the matrix 


d 

A{u,uj) := ^^u}kA^{u) ( 1 - 3 ) 

k=l 

has n real eigenvalues Ai(m, w),..., An(u, w) and a complete set of left (resp. right) eigen¬ 
vectors li{u,u}) = .. ,lin{u,uj)) (resp. ri{u,uj) = (rii(u,w),...,r„i(u,w))^) (i = 

1,..., n). It is assumed that Xi{u,uj) {i = 1,..., n) are smooth with respect to u and lo, while 
li{u,io), ri{u,io) {i = 1,... ,n) are all smooth functions of u for any given a; G One may 

normalize the eigenvectors so that 


k{u,uj)ri^{u,u}) = 5iii, = 1,... ,n. 


(1.4) 


It is well known that for any given suitably smooth initial data uq, the system (HH) has a unique 
local smooth solution mms], but in general the solution develops singularities in hnite time 
even for small and smooth initial data [ana US]. 
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It has been known that there are two structures that can prevent the finite-time singularity 
formation. One of them is the (weak) linear degeneracy for the homogeneous systems, which 
was firstly introduced in ID study (see [3] and m)- The other one is the dissipation effect 
induced by inhomogeneous terms. A well-known assumption in this aspect is given by the strict 
dissipation, which requires a damping term to enter into each of the equations of the system, 
see m- Notice that for some physical problems, this condition is too strong, while in many 
applications and as the interest in this paper Q{u) takes the form 

Q{u) = (0,..., 0, Qr+i{u ),..., Qn{u)f', (1.5) 

which seems to produces the dissipation effect only in the last n — r equations. In this case, the 
system (HU) is referred to as a partially dissipative hyperbolic system in many literatures, such 

as [I1I21IIIE5]. 

A constant vector u* £ M”' is called an equilibrium state of the system (11.11) if 5(u*) = 0, 
or equivalently, (5 (m*) = 0. The main subject of this paper is to search for suitable structural 
conditions, under which (HU has a unique global smooth solution near u^- Without loss of 
generality, one may take u* = 0, he., (5(0) = 0. One may also suppose that G(0) = 0 and 
V„G(0) = In- A useful assumption (see [25]) is the following entropy dissipation condition: 

(211) The matrix ^f^(O)^ is invertible. 

( 2 1 2 ) There exist a strictly convex smooth entropy function f)(G) and d smooth entropy flux 

functions {k = 1,... ,d) for the system (jl.ll) in terms of G, such that for all G, 

{VGi^HG)f = {VGfi{G)fVGiF\u{G))). 

This implies that if one defines r]{u) = fi{G{u)) — f)(0) — Vgv{ 0) ■ G{u) and 'ij:^{u) = 
■ii)^{G{u)) — Vg??(0) • F^{u), then r/(0) = 0, Vur]{0) = 0, V^r/(0) is positive definite and 

= {VuV{u)fA\u), k = l,...,d. (1.6) 

( 2 1 3 ) There exists a positive constant Ce > 0 such that for all u, 

V„r/(M) • Q{u) < -Ce\Q{u)\'^. (1.7) 

It is clear that the entropy dissipation condition (2ti)-(2l3) is too weak to prevent the formation 
of singularities, for instance, a decoupled system composed of an inviscid Burgers’ equation and 
a damped linear equation satishes ( 2 li)-( 2 l 3 ), but generally would develop singularities in a 
finite time for non-degenerate initial data. 

A natural supplementary coupling condition, the Shizuta-Kawashima stability condition, or 
“Kawashima condition” for short which was first formulated in na |23| for the hyperbolic- 
parabolic systems with a source term, was introduced in mm- This coupling condition im¬ 
plies that due to the wave propagation, the partial dissipation term indeed produces a kind of 
complete dissipation effect with respect to all the components of the solution, which yields the 
global existence of the small smooth solutions. The Kawashima condition has several equivalent 
formulations, and one of them reads as 

(.h) The kernel of the Jacobian Vu(5(0) contains no eigenvector of A(0,a;) for any lo G 

Under the assumptions ( 2 ti)-( 2 t 3 ) and (.^), [Z1 [25| prove the global smooth solutions for small 
initial data. See also [n El [n [m ESI for more discussions and |5| for the result on global 
entropy solutions to systems satisfying Kawashima condition with zero mass initial data. 

It should be pointed out that although the name partial dissipation is used in this case, all the 
components of the solution are dissipative due to the Kawashima coupling condition (R). On 
the other hand, any other kind of partial dissipations that certain part of the solution may not 
be linearly affected by the dissipation and hence is non-dissipative, which can be represented by 
the violation of the Kawashima condition (.^), is an interesting structure to study (see remarks 
in m Eoi 125] 1. In fact, there are many physical models that the Kawashima condition (.^) 
is violated, and the eigen-family that violates it is linearly degenerate, such as the damping 
full Euler system of adiabatic flow (see, for instance, ( 8 ] Elj) and the thermal non-equilibrium 
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flow (see [Ml 1271 [28]). Some works have discussed several aspects on the system violating the 
condition (.^), such as [M] for special structures and [18] for partial strict dissipations in ID. 

In this paper we will continue the study on the system with d > 2 and r >2 (and hence 
n > 3) in (11.51) violating (K). More precisely, we assume 

( 2 I 4 ) For each cj G one and only one family of right eigenvectors of ^( 0 , 0 ;) are not 

contained in the kernel of the Jacobian V„(5(0), that is, possibly after relabeling, 

VuQ{0)rj{0,uj) ^ 0, j = 2, ...,n, (1.8) 

but 

V„Q(0)ri(0,u;) =0. (1.9) 

One may refer to the assumption ( 2 I 4 ) that the last n—1 right eigenvectors satisfy the Kawashima 
condition and the first right eigenvector violates it. Note that ( 2 ti)-( 2 t 4 ) is still too weak to 
prevent the formation of singularities, and some additional structural conditions are needed in 
order to guarantee the global existence of the small smooth solutions. The goal of this paper is 
to identify one set of such conditions and prove the corresponding results. 

For notational simplicity, and are used to denote the usual and Sobolev spaces on 
M'^. The notation is defined by A^/ := (|‘^|^/)^, where " is the Fourier transform and (•)^ 
is its inverse. C denotes the positive constants, whose value may change from line to line, and 
f <g IS used for / < Cg. 


2. Main results 


As mentioned in Section [U the aim of this paper is to impose certain structural conditions 
in supplement with the assumptions (2ti)-(2l4) so that the system (II.ID admits a unique global 
smooth solution for the small initial data. First, for the technical point and as suggested by 
many physical applications, one may assume 

(*B) The first family of eigenvectors satisfy an isotropic condition: 

li{u,co) = li{u) and ri{u,Lj) = ri{u), Vo; G (2-1) 

In the spirit of the normalized coordinates method developed in [IlKIT] for the one dimensional 
quasilinear hyperbolic systems, thanks to the isotropy of ri in (12.ip . the partially normalized 
coordinates can be introduced as follows. After the linear transformation performed in Appendix 
S let u*^) be the 1st characteristic trajectory passing through the point u^, which is 

defined by 

I^ n(ix^^^(s,w°)), ^2 2) 

u^) = vP. 

For each point u = (fti,..., Un)^ , map it to u = (ui,..., Un)^ by 

u = u{u) := uP^{ui, {0,U2, ■ ■ ■ ,Un)'^)- (2.3) 

It will be verified in Section [3] that the mapping (12.3D is a local diffeomorphism near ft = 0. This 
then allows to define the partially normalized coordinates u by the inverse mapping u = u{u) 
for each point u near tt(0) = u^^^(0,0) = 0 in the phase space, and hence the system (II.ip can 
be equivalently reformulated in terms of u. This reformulated system will be studied in details 
in Section (3] Throughout the paper the following notational convention is applied for vectors: 
for a vector / = (/i,..., fn)'^ G M”, denote 

= (/2,..., fnf e f = (/2,..., frf G and f = ifr+1, ...Jnfe (2.4) 


As one shall see, the decomposition ()2.4p is closely related to the structure of the system (II.ip 
under assumptions ( 2 li)-( 2 t 4 ) and (IB). 

Since the assumption (Stq) implies that the first right eigenvector violates the Kawashima 
condition, the natural way is to supplement certain degeneracy conditions with respect to the 
first eigen-family. One may use the following weak degeneracy conditions: 
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(2nDi) Q{u) is weakly degenerate along the first characteristic trajectory that passes throngh 
0 : 


(5(n^^^(s,0)) = 0, V |s| small, 


which is equivalent to 


(V„Q(n)ri(u))|„^„(i)(^o) = 0, V|s| small. (2.5) 

{WDi) The first eigen-family is weakly linearly degenerate (see [161117] l: 

Ai(n^^^(s,0),a;) = Ai(0,w), V |s| small, 
which is equivalent to 

(V„Ai(u,w) • ri(n))|„^„(i)(^ Q) = 0, V |s| small, Vw e (2.6) 

It is easy to check that the conditions (2IhDi)-(2IT1D2) are necessary since if any one of them 
is missing, the solution generally will blow up in finite time. The typical examples can be 
constructed in the spirit of the Riccati equation and the Burgers’ equation (see |16] and |18] for 
instance). 

The first result of this paper is that under the assumptions (2li)~(2l4), (53) and (21hDi)~ 
(WDz), the system (II.ip admits a unique global solution for each suitably smooth and small 
initial data. 


Theorem 1. Suppose that the assumptions (2ti)-(2l4), (53) and (2n(i)i)-(2ItD2) hold. Let 
d > 3, i > d/2 + 1 and 1 < p < 2, 1 < p < d/2. Denote p* = m.m{2, dp/{d — p)} and 
sj = d{l - 1/p) + 1. 

There exists a sufficiently small eo > 0 such that if 

Kq ■= ll^oll//^ + II^oIIlp II^oIIlp* — ^0) (2-7) 

then there exists a unique global solution u(t, x) to the system (jl.ip satisfying 


+ Ah^(r) dr < 

(2.8) 

Moreover, for all s G [0, s^] n [0, d/2), 


\\Affi^{t)\\^,_^<CKo{l + t) ^(p 2 ) 2 

(2.9) 

and for all s G [0, — 1] n [0, d/2 — 1), 


\\Affi°{t)\\j^,_^<CKo{l+t) 2 ) 2 2 , 

(2.10) 


Remark 2.1. Note that if d > 5, then one can take p = 2 in Theorem{l\ which yields the unique 
global solution to the system a for the small initial data in with i > d/2 + 1; but if d = 3 
or 4, in order to obtain the global solutions one needs the low frequency assumption of the initial 
data as stated in (1221) with 1 < p < d/2 (if p > 2d/{d -|- 2) and hence p* = 2, then one needs 
only the low frequency assumption of Uq). 

To improve the result of Theorem[Tl one may use the following stronger degeneracy conditions: 
((Di) Q{u) is degenerate along each first characteristic trajectory: 

NuQ{u)ri{u) = 0, Vu. (2T1) 

((D 2 ) The first eigen-family is linearly degenerate (see |15jl: 

VuXi{u,oj) ■ ri{u) = 0, Vtt, Vw e (2-12) 

Then the second result of this paper can be stated as follows. 
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Theorem 2. Suppose that the assumptions (2li)-(2l4), (55) and (T>i)-(Ti 2 ) hold. Let d > 3, 
I > d/2 + 1 and 1 < p < 2, require further that p < 2d/(2(3 — i) + d) if i < 3. Denote 
^2 = min{d(l — 1/p) + 1,£ — 1}. 

There exists a sufficiently small eo > 0 such that if Kq < eq, then there exists a unique global 
solution u{t, x) to the system (II.ip satisfying 

\Ht)\\Hi + J^ (^\\u°{T)\\]^t + \\^u^{r)\\%-i^dT <CKl (2.13) 

Moreover, for all s G [0, s^] [Oj d/2 + 1), 

||A^u"(t)||^,_, < CKo{l + (2.14) 

and for all s G [0, — 1] n [0, d/2), 

||A^il°(t)||^,_, < CKoil+t)-^ir-^)-^^-\ (2.15) 

Remark 2.2. Note that now one can take p = 2 in Theorem\^for any d >3 which yields the 
unique global solution to the system dm) for the small initial data in with i > max{d/2+l, 3} 
and without any low frequency assumption of the initial data. 

Note that the result of Theorem [2] excludes the case d = 2. To obtain the global solution to 
the system dnD in two dimensions, one more assumption should be made on the structure of 
the system: 

((D 3 ) The first eigen-family gets no linear effect from the source term: 

^i(0)V„Q(0) = 0. (2.16) 

And the third result of this paper is stated as follows. 

Theorem 3. Suppose that the assumptions ( 2 ti)-( 2 t 4 ), (55), ((Di)-((D 2 ) and {T)f) hold. Let 
d>2,i> d/2 + 1 and 1 < p < q < 2, 1 < p < q < d, require further that p < 2d/ (2(3 — ^) -|- d) 
if £ < 3. Denote S 3 = min{d(l/2 -|- 1/q — 1/p) + l,d(l — 1/p) + 2,£ — 1}. 

There exists a sufficiently small Eq > 0 such that if 

Kq ■= II'^oIIh^ + II'^oILp + II'^oILp* + ll^o,i|lx,9 < £ 0 , (2-17) 

then there exists a unique global solution u(t, x) to the system (jl.ip satisfying 

\\uit)\\Hi + (||'w°( 7 ')|Ih^ + dr < (7.^^. (2.18) 

Moreover, for all s G [0, S 3 ] H [0, d/2 -|- 1), 

||A^u"(t)||^,_, < CKoil + ( 2 . 19 ) 

and for all s G [0, S 3 — 1] n [0, d/2), 

||AV(t)||^,_, < CKo{l+t)-^ii^-^)-^^-K ( 2 . 20 ) 

Remark 2.3. It is verified in Section\JU that the damping full Euler system of adiabatic flow sat¬ 
isfies all the assumptions ( 2 li)-( 2 t 4 ), (55) and ((Di)-((i) 3 ) and hence Theorem\^can be applied. 
This improves the result of [24] in 3D and produces the first result of global unique solution in 
2D which seems necessary to require the low frequency assumption of the initial data. However, 
our results could not cover the multidimensional thermal non-equilibrium flow |28| . 

Now let us illustrate the main difficulties encountered in proving these results of Theorems 
HHSl and explain our strategies to overcome them. Since the components of the solution behave 
totally distinctively as already seen in the theorems, it is essential to divide the system, at least 
in the linear level, into the dissipative part and the non-dissipative one. Different from the 
specific models which usually have certain physical variables to give reasonable ways of doing 
this, for the general system considered in this paper these two parts are generally strongly 
coupled in an implicit manner. In [18], the methods of normalized coordinates and formulas 
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of wave decomposition are used to identify the two parts for the one dimensional systems with 
partial strict dissipation. However, for our multidimensional systems, these methods, which 
are originally designed only for one dimensional problems, are just not applicable. Meanwhile, 
when discussing the dissipation structures given by the entropy dissipation condition together 
with Kawashima condition instead of the strict dissipation, the way of applying these methods 
would generally destroy the structure to identify the faster decaying variables from the slower 
ones. To get around this difficulty, the partially normalized coordinates (12.31) is introduced 
under the assumption (55). It is shown in Section [3] that by reformulating the system (II. ip in 
terms of this new coordinates into the system (13.81) . the two parts are successfully identified, 
z.e., the dissipative part ()3.26p for and the rest part (|3.25p for the non-dissipative 

ui. Also, those structural conditions can be reformulated in terms of new variables in a more 
applicable way. 

The main part of proving Theorems [TH3] is to derive the energy estimates. In the following, 
we explain the strategy by first dealing with the case under the assumptions of Theorem [U 
Theorems [2H3] are proved in the same fashion and we only need to explain the additional 
arguments for the proof. It is crucial to analyze the structures of the linearized system and the 
nonlinear terms. It is with the help of the partially normalized coordinates that the linear part of 
(I3.26P is shown to be dissipative in the sense of satisfying the entropy dissipation condition and 
the Kawashima condition. By putting aside the nonlinear effect, one can follow the arguments 
of [TIES] to derive the basic estimate by using the entropy dissipation condition ( 2 li)-( 2 l 3 ) 
as well as the high order energy estimates for v!' which include the damping dissipation of 12 °; 
the partial Kawashima condition ( 2 I 4 ) can be applied to recover the degenerate dissipation of 
ft*". The remaining difficulty is to handle with the nonlinear terms. It is again with the help 
of the partially normalized coordinates that under the degenerate assumption ( 2 B(Di) there is 
no nonlinear term resulting from the source term composed solely by the non-dissipative part, 
that is, those nonlinear terms are of form 0{l)uu^, at least. This feature provides relatively 
nice structure so that we can complete the estimate of the dissipative part . 

The next natural step is to derive the energy estimates for the remaining non-dissipative 
part ui- Under the weakly linearly degenerate assumption (SIhDz), there is no nonlinear term 
resulting from the convection term in ()3.25l) composed solely by the non-dissipative part. How¬ 
ever, there would be a loss of derivative if one directly performs the energy estimates by using 
(I3.25P . To get around this difficulty, we return back to the original full system (|3.8p and resort 
to generalize the method of wave decomposition (see M) to our multidimensional case. The 
isotropy assumption (55) will be used again. After the delicate analysis of the nonlinear term 
under the degenerate conditions ( 2 IfDi)-( 2 IfD 2 ), we can finish the full energy estimates as 


\u 


Hi 


f 


+ ■ ■ ■ dr. 


( 2 . 21 ) 


Note carefully that we cannot control the time integral in (12.211) by the dissipation estimates. 
To bound this time integral, we will turn to derive a strong enough decay rate of 
This step is extremely delicate. We will need to highly explore a refined linear decay estimates 
which requires weaker low frequency assumption of Uq than Uq, the stronger degeneracy of Q'" 
than and the faster decay of than u*". The technical point thus lies in the way of balancing 
these despairs. Taking into account these facts, by a bootstrap estimates, which implies that 
the decay of higher order norms can be deduced from the decay of lower order norms, and a 
use of the Duhamel formula, which implies that the decay of lower order norms is bounded 
in terms of the decay of higher order norms, we are able to derive a decay of the solution as 
stated in Theorem [T] in a recursive way by the smallness of the solution. It then follows from 
the Sobolev interpolation that is integrable in time as long as the parameters p and d 

satisfy the required restriction given in the theorem. We remark here that the introduction of 
the fractional derivative in the study enlarges the available range of p and d. This closes the 
energy estimates and then completes the proof of Theorem [TJ 
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To prove Theorem[2l it is observed that under the stronger degenerate conditions (!l)i)-(S 2 ), 
the nonlinear terms behave better. For example, under (Si), nonlinear terms resulting from the 
source term are of the form 0{l)uiP, or so. Based on the stronger degeneracy, we can improve 
the full energy estimates as 

^ + ---dT. (2.22) 

Since iP decays at a faster 1/2 rate than if, this ultimately enlarges the range of p and d as 
stated in Theorem [2j Note that we can mostly deduce from (|2.15p that decays at the 

rate of (1 + fQj. g^j^y g > q, which restricts Theorem [2] only hold for d > 3. 

The main goal of Theorem [3] is to include the case d = 2. Note that even in view of the 
linear decay estimates decays at the rate of (1 + and then we need to 

require the low frequency assumption of when d = 2. However, due to the nonlinear effect 
we cannot achieve this rate, and Theorem [2] basically tells that would decay at the 

same rate of (1 + £qj. p ££ £g further additional assumption. Thereby, 

the assumption of ho,! comes into the role. The idea is to derive the estimate of ui 
so that we can improve a bit some nonlinear estimates when using the Duhamel formula, for 
example the estimates of However, we can not derive the L'^ estimate of ui by using 

(I3.25[) , and we again need to return back to the original full system ()3.8p and introduce the wave 
decomposition to consider instead the first wave vi. It is the derivation of the estimate of vi 
that we need the assumption (03). With the help of the estimate of hi, we can improve the 
nonlinear estimates so that we can derive the decay rates of the solution as stated in Theorem 
[3l Note that now the decay rate of is improved to be (1 + which allows 

the case d = 2 for q < 2. One may then notice the necessity of the assumption of uo,i. 

The rest of the paper is organized as follows. In Section [3l partially normalized coordinates 
will be used to reformulate the system (tni) and those structural conditions. In Section H] 
the structure of the reformulated system is analyzed carefully as the foundation for the further 
analysis. The delicate energy estimates is performed in Sections [SHSl the low frequency estimates 
for two dimensional case is presented in Section [71 and the refined decay estimates is given in 
Section [8l Finally, those estimates are combined in Section [9] to prove Theorems [TH3l 

3. Partially normalized coordinates 

After a linear transformation if necessary (see Appendix [A] for details), one may assume 

^(0)=0, g= l,...,r; z = l,...,n, (3.1) 

and 

n(0) = ei, (3.2) 

where (z = 1,..., n) stands for the ith unit vector in M”. 


The first step is to verify that the nonlinear mapping u = u{u) in the phase space, dehned 
by (|2.3I) and (|2.2p . is invertible near u = 0. Denote the Jacobian matrix of the transformation 
u = u{u) by 



JW := ||(u). 

(3.3) 

By the definition, it holds 



II 

II 

(3.4) 

and 

du , . , , ., 

(3.5) 


— (h) = ri(u(h)). 

Hence, 


(3.6) 

This together with ()3.2p yields 

J(0) = (ri( 0 ),e 2 , ...,en) = In, 

(3.7) 
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which implies that the transformation u = u{u) is a local diffeomorphism near u = 0. In the 
partially normalized coordinates u of the phase space, the system (11.11) (and (11.21) 1 is transformed 
into 

' d 

dtu + ^ A^{u)uxt, = 
k=l 

1,11 |t=0= Uo- 

Here the coefficient matrices 

A^(u) = J~^(u)A(u(u))J(u), k = l,...,d (3-9) 


(3.8) 


and the inhomogeneous term 

Q(u) = J~^(u)Q(u(u)). (3.10) 

Now it is important to study the structure of A^{u), which will exhibit the advantage of our 
partially normalized coordinates u. First of all, it is direct to check that A^{u) possesses the 
following expression 

A^{u) = A{u,ek), fc = l,...,d (3-11) 

with Gfc the /cth unit vector in and 

A{u,oj) := J~^{u)A{u{u),uj)J{u), io € (3-12) 

where the matrix A{u,uj) is defined by (11.31) . Therefore, it suffices to study the structure of 
A{u,uj). Recalling that li{u,uj) (resp. ri{u,uj)) is the ith left (resp. right) eigenvector of A{u,uj) 


with the corresponding eigenvalue Xi{u,uj), i = 1,... ,re, one could then choose 

Xi{u,uj) := Xi{u{u),uj), (3.13) 

ri{u,uj) = {fii{u),... ,fni{u))^ := J~^{u)ri{u{u),uj), (3.14) 

li{u,uj) = ... ,lin{u)^ := li{u{u),uj)J{u) (3.15) 

as the eigenvalues and a complete set of left (resp. right) eigenvectors of A{u,uj), that is, 

li{u,uj)A{u,Lj) = Xi{u,u})li{u,Lj), i = (3.16) 

A{u,uj)ri{u,uj) = Xi{u,uj)fi{u,u}), i = (3-17) 

From the normalization (m, it follows 

li{u,uj)rif{u,uj) = i,i'= 1,... ,n. (3.18) 

Moreover, by the assumption (*B) and (|3.5p . respectively, for any ui S 

li{u,uj) = li{u{u))J{u) = li{u) (3.19) 

and 

fi{u,uj) = J“^('u)ri(tt(?2)) = ei. (3.20) 

Combining (|3.18p - (l3.2Up . it follows 

[ 11 ( 11 ) = 1 and (11,0;) = 0, j = 2,..., n, Vo; € (3-21) 

More importantly, taking 1 = 1 in (13.171) and noting (13.201) . leads to 

Hii (H, o;) = Ai (H, o;) and Hji (H, 0 ;) = 0, j = 2,..., n, V 0 ; £ . (3.22) 

Throughout the paper the following notational convention are used for matrices: for an re x re 
matrix M = denote 

M« = (Afi 2 ,..., Mi„), = (M,y); .,^2 and . (3.23) 


Then by (I3.22p . one has the following structure of A(u,uj): 

Till o;l - I“(ll,w) 

0 A^{u,oo) 


(3.24) 
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Accordingly, the system (|3.8p indeed takes the form 

d d 


dtui + '^Xi{u,ek){ui)xk + = Qi{u), 


k=l 

d 


k=l 


dtu'’ + ^A’^’\u)ui^=Q\u) 


(3.25) 


(3.26) 


k=l 


Note that the system (13.81) is split into two subsystems in some sense. We will see that the 
subsystem (j3.26l) of is dissipative in the sense of satisfying the entropy dissipation condition 
and the Kawashima condition and the main part of the subsystem (I3.25P of the non-dissipative 
ui is indeed (weakly) linearly degenerate. This structure property is crucial for the whole 
analysis in this paper. 

Recall that from the assumption (*B), the isotropy (|3.19p - (l3.20p hold in the partially normal¬ 
ized coordinates, and we restate them as follows 

(*B) The first family of eigenvectors satisfy an isotropic condition: 

li{u,u}) = Ii{u) and fi{u,oj) = ei, (3.27) 

The next task is to derive the corresponding forms of the assumptions (2ti)-(2t4), (SlhDi)- 
(WDi), {Di)-{Di) and (T>3). First, recalling (I3.10p and (11.51) it trivially holds that (5(0) = 0 
and 


'^Jqi{u)Qi{u) = 0, q = l,...,r. 


1=1 


Moreover, by (13.101) and (13.71) . one has 


diii 


( 0 ) = 


dJ 


-1 


dui 


-(u) 


u=0 


Q(o)+j-'{o) (x: ##(«(»)) 


\i'=i 


dui duj! 


u=0 


- £(»)■ 


(3.28) 


(3.29) 


Hence the assumption (2ti) is translated to be 
(2ti) The matrix ^Vu(5^ (0) is invertible. 
Next, set 

then 


fj{u) = viuiu)), 

f)(0) = 0 and Vui?(0) = 0. 


Moreover, since 
d^fj 


_(i) = + x: (- 1 ^( 0 ) 

' dUrdUrl dUi dUi! V dUr- dUr- 


dur , _, du. 


dmdur 


r,r'=l 


r=l 


drj 
dur 


drj 
dur 


d'^Ur 

dujduj 


(3.30) 

(3.31) 


(3.32) 

(3.33) 


by (13.71) it holds 

V|i?(0) = J^(0)V,%(0)J(0) = V57(0). 

Set '0^(u) = ip^{u{u)), then the assumption (Stz) leads to 

(Stz) The strictly convex smooth entropy function fj^u) and the d smooth entropy flux func¬ 
tions 'il:^{u) {k = 1,... ,d) satisfy that for all u, 

= {VuV{u)fA^{u). (3.34) 

Now, noting (13.71) . (13.101) . (I3.14p and (I3.29p . it follows from the assumptions (Sis) and ( 2 I 4 ) 
that 

( 213 ) There exists a positive constant Cg > 0 such that for all u, 

^ud{u) ■ Q{u) < -561^(11)1^. 


(3.35) 
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( 2 I 4 ) The last n—1 right eigenvectors of ^(0, u) are not contained in the kernel of the Jacobian 
VnQ(O) for any uj £ that is, 

VuQ{0)fj{0,uj) 0, j = 2, (3.36) 

but 

v^g(o)fi(o) = 0. 

For the assumptions ( 2 iriDi)-( 2 iriD 2 ), we notice that the dehnition (12.31) directly leads to 

tt(uiei) = 0). (3.37) 

Then by (|3.10l) . (|3.13l) and (I3.27p . one can express them as 

( 21 hDi) Q{u) is degenerate along the hrst characteristic trajectory passing through 0 : 

—^z^('Uiei) = 0, Vluil small. (3.38) 

OIL I 

{WDi) The first eigen-family is weakly linearly degenerate: 

7 ^(niei, w) = 0, V luil small, Vw £ (3.39) 

OUi 

Similarly, the assumptions (T>i)-(T> 2 ) are translated to be 
((Di) Q{u) is degenerate along each first characteristic trajectory: 

^^^(u) = 0, Vu. (3.40) 

(Til) The first eigen-family is linearly degenerate: 

r)X 

-—^(u,u}) = 0, Vu, Vcj £ (3-41) 

OUi 

At last, by (13.7p . (I3.19p and (I3.29p . the assumption ((D 3 ) can be expressed as 
((D 3 ) The first eigen-family gets no linear effect from the source term 

/'i(0)VsQ(0) = 0. (3.42) 

Remark 3.1. Since the assumptions (2ti)-(2t4) and (55) are made in all the three Theorems 
[2H5l (2li)-(2t4) and (55) are always supposed to be hold in the rest of the paper. 

4. Structural analysis 

In this section, useful algebraic consequences are established for the structural assumptions 
( 2 ti)-(a 4 ), (i^i) and (3i). 

4.1. Structure of the source term Q. In this subsection the source term Q{u) is analyzed 
under the assumption (2IlDi) or ((Di) by the Taylor expansion in an appropriate way. 

4.1.1. The expansion ofQ{u) under (QTTii).' First, assume (5IlDi), then it holds that Qiuiei) = 
J“^('Uiei) J(0)(5(0) = 0. Using this fact, one may first expand Q{u) around uiei, that is, hxing 
ui and viewing Q{u) as a function of = {u 2 , ■ ■ ■, Un)'^ and then expanding it around v}’ = 0 , 
to have 

n QQ n 

Qi{u) = -^{uiei)ujf + 

f=2 j',j"=2 


i = 1,... ,n. 


(4.1) 
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Expanding further the functions of ui inside the first summation in (|4.ip around ui = 0, one 
can obtain 


j'=2 ^ j'=2 ^ ^ j\j''=2 


(4.2) 


j'=2 


n. 


Now examine the hrst two summations in (|4.2p . First, it follows from (13.291) . (13.11) and (II.5p 
that 


dQi dQq 


dUn 


(0) = -^(0) = 0, g = 1,..., r; i = 1,..., n. 


(4.3) 


On the other hand, by (j3.10l) . it holds 

OUi UUi> OUi OUi> OUi'C/Ui 

Hence, by the fact Q(0) = 0, (|4.3p . (13.71) and (II.5p . one has 


(4.4) 


d^Q 




dUqidUq” 

Denote the constant matrix 


and for i, i', i” = 1,... ,n, 


^ {o) = —^{Qqiuim 


dUqidUqii '' 


= 0, q,q',q' = 1,... r. 


u=0 


0 = ViiQ(O) 


T rni 

OUiiUUiii 

Then (|4.3I) and (|4.5I) imply that 

Qqi = &iq = 0, z = 1,... , n; g = 1,..., r, 

'^qq'q"=0, q, q', q" = 1, ■ ■ ■ T. 

Accordingly, it follows from ()4.2I) that 


(4.5) 


(4.6) 


(4.7) 


(4.8) 

(4.9) 


Qqi^) = X] '^qp'iUp'Ui+ ^ Yqjij//{u)ujiUjfi+ '^ Zqjfii{u)ujiul, q = l,...,r (4.10) 


p'=r-\-l 


r,j"=2 


j'=2 


and 


Qp('i^) — ^ ^ ©pp'flp' T ^ ^ Ypj/p/(^ujUj/Ui//, p — r + 1,..., n, 

p^=r+l j'=2i"=l 

where Ypjii{u) = Ypiji + Zpjiii{u)ui. 

Remark 4.1. By the matrix 0 defined in (|4.6p . (2ti) is equivalent to 

det 0° 0. 


(4.11) 


(4.12) 
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4.1.2. The expansion of Q{u) under (2)i). Since (iDi) implies (QUSi), the expansions ()4.10p - 
(14.lip and all the computations in the previous subsubsection hold in this case. But under the 
stronger assumption (iDi), it is reasonable to expect a better expansion. 

Again, using the fact Q(uiei) = 0, one may first expand Q{u) around uiei to have 


Qi{u) = 


n 


dQi 

^ duji 
j'=2 ^ 


{uiei)uj> + ^ 


1 d‘^Qi 


2 duiiduii 
j',j"=2 ^ ^ 


-{uiei)ujiUjii 


n 

i'll" j*=2 


i = 1,... ,n. 


(4.13) 


Due to the assumption (iDi) and the fact Q(uiei) = Q{0) = 0, it holds 

Q{uiei + UgCg) - Q{uiei) =J~^{uiei + Uqeq)J{Uqeq)Q{Uqeq) 

— J (filCl + UqCq') J (Uqeq') i^QiUqOq^ Q(0)^ . 

This together with (14. 3 p yields 


^(uiei) = J ^(£(iei)J(0)^(0) = 0, g = l,...,r. 

OUq aUq 


Thus, from (I4.13h . one can get 


3.W = E + E 


UjiUjii 


p'=r+l 


j'd"=2 


i':i" J*=2 


(4.14) 


(4.15) 


(4.16) 


Expanding further the functions of ui inside the first two summations in ()4.16l) around ui = 0, 
one has 


dQ 


Qi{u) — ^2 Qu 


,,l[U)Up>Ul 


p'=r-\-l 
1 


p'=r+l 


Y 




(4.17) 


-(0)UjrUj» + '22 Y^ {u)UjiUji/Ui* , Z = 1, . . . , U. 


2 dupdup 

j/J"=2 ^ J jij"=2i* = l 

Recall 0 from (j4.6p and T from (14.7p . and note ()4.8p and ()4.9li . then it holds 


Qq(.u) = '22 X] + Y X] 9 = 1,---T (4.18) 

p'=r-\-li"=l j2j"=2i* = l 


and 

n n n 

Qpi^) = Y ®pp'V+ Y Y2'i{u)up>ui+ '22 P = r + l,...,n, (4.19) 

p'=r-\-l p'=r-\-l j2j"=2 


where Y' ,.„{u) = \Yqp,jn and Y'.,.„{u) = ^Tppp, + Z'{u)u 


4.2. Structure of the subsystem (|3.26l) . In this subsection the subsystem (I3.26P is analyzed 
under the assumptions ( 2 ti)-( 2 t 4 ). 
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4.2.1. Symmetrizer. The first task is to search for a positive definite (n — 1) x (n— 1) symmetric 
matrix to symmetrize the subsystem (I3.26p . Indeed, 


A^’\u) : = 






dujduj 




dG 


-1 


dujduji 


(u) 


(4.20) 




is just the desired symmetrizer. 

Lemma 4.2. is symmetric and positive definite near u = 0. Moreover, 

Proof. It is direct to check that A^’^{u) is symmetric. Meanwhile, by (jd.dll) . it holds 

I°>(0) = (V|r)(0))^ 


(4.21) 


(4.22) 


On the other hand, according to the assumption (Stz), fj is strictly convex, then it follows from 
(IT^ that M°’‘'(u) is positive definite near u = 0. Moreover, it follows from (13.Op and (I3.34p 
that, setting F^{u) = F^{u{u)), 


0 = 


d ( 


d ( 


dui' \ dui I 

d^f] 


a 

rj=l 


dujtdur V du 


dui \ dui! I 

dG\~^dFl^ 


+ 


df] d f f dGY^\ dFj;^ dp f dG 


rl 


duj dur dujt 


du J 


rl 


dui dur V du 


-1 Q2pk 


-E 

r,l=l 
n / 


d'^p (dG\~^dFj^ ^ dp d 


duidur J rl dUi 


du ) 


-F 


rl 


+ 


duidui! ^ 
-1 o2 pfc 


dp fdGY 




dur dur \ du J ,.1 durdF 


d^p ^ dp d 


^ i dUiidUr \du ) , dUr'^^'^ dUr dUi! 

,r',1=1 \ \ / iL 


du ) 




rl 


dUr' 


d^ ^ dp d ((dG\ 


\^duidur \du ) dUr' dUr dUi ^ ^ 


: 


dUr' 


I r^o ~ n 

^ duudur ^ 

. r=l * ^ r,r',l=l 


Y. - 


df, /day' ff^G, 

dUr \du ) dUildUr' 


dp fdGY d'^Gi 


, r=l 


dUfdUr^ 


^ dUr I j ; dUidUr’ 


This together with (I3.22p gives rise to (I4.2ip . 


□ 


4.2.2. Dissipative structure. In this part the linear dissipative structure is investigated for the 
subsystem (I3.26p . which yields the damping effect on the component uP of the system p3.8p . 


Lemma 4.3. Recall A'^’^(O) from (I4.20p and 0 from (|4.6I) . Then it holds 


Al°’‘'(0) 


0 0 
0 0° 


+ 


(e”)^) '^"•‘(o) = - (I 


0 ^ 
0 M° 


for a constant [n — r) x (n — r) symmetric matrix M°, which is positive definite 


> Cmln-r, 


(4.23) 


(4.24) 


for a constant Cm > 0. 
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Q 0 D )) block accordingly V„r/(0) = 

V,'r 7 ( 0 )V„Q( 0 ) + (VuQmfV^m = 


Proof. We employ a modification of the proof of Theorem 2.1 of [25]. Recall that V„(5(0) = 

0 0 \ ..' 

i ?2 -R 4 


. It follows that 


+ 




Bi Bi\ fO 0 
B2 R4 ) VO 0° 

0 

( 0 °)^R 2 R40° + (0°)^54 


0 0 
0 ( 0 °)^ 


Bi 

B 2 Bi 




(4.25) 


Since {yuV{u))'^Q{u) achieves its local maximum at h = 0 by ( 2 I 3 ), one gets that V|? 7 ( 0 )Vji( 5 ( 0 )+ 
(^uQ(0))^V?fj(0) is semi-negative dehnite. It then follows from (I4.25P that 


V,'i 7 ( 0 )V„Q( 0 ) + (VuQ(O)fV^m = - 


0 Jl 

0 M° 


(4.26) 


where = —i? 40 ° — (Q^j'^Bi is symmetric and semi-positive definite. 

Note that (14.231) follows by (14.261) and (I4.22p . It remains to prove that M° is indeed positive 
definite. Note that it is equivalent to prove that is invertible. Suppose that has a zero 

eigenvalue, one may take a corresponding right eigenvector it° / 0 and set it* = ( .Pn V Then 


by a Taylor expansion, since Q(0) = 0 = Vu^(O), 

0 >yQy) + Ce|Q(it*)p 

(V|?7(0)ViiQ(0) -h (ViiQ(0))^V|f7(0))it* -h Ce\VuQ{0)uy -h o(|it*|^) 

= - (it*)^2M°it* -h Ce|0°it*p o(|it*p) 

=Ce\eyy + o{\uy), 

which yields the contradiction since 0° is invertible by (2li), if it* has been taken sufficiently 
small. □ 


4.2.3. Partial Kawashima condition. Finally, the partial Kawashima condition inside the as¬ 
sumption ( 2 I 4 ) for the subsystem (I3.26p is explored, which implies that the component it*" of 
the system (|3.8p is indeed affected by the dissipation linearly. 

Lemma 4.4. For any oj G recall A^(0,a;) from (I3.24p . Then there exists a constant > 0 
and an (n — 1) x (n — 1) skew-symmetric real matrix K = K{ijj) G C'°°(§'^“^) with 

K{-u) = -K{uj) (4.27) 

such that 

K{uj)A'’{0,u;) - (^A\0,u;)Y K{uj) > 2cfc/n-i - 2diag(0, R_i), V w G (4.28) 

Proof. Consider the following linear system 

I 0 ,b( 0 )^^b rV (4.29) 

k=l ^ ^ 

Due to (|4.3I) and (|3.36l) . it holds that f^(0,a;) 7 ^ 0 and hence 

(0 i = 2,...,n. (4.30) 


On the other hand, by (13.171) and (13.241) . 



is a set of right eigenvectors of the 


matrix A^{0,u}) with Xj{0,uj) (j = 2 ,...,n) as the corresponding eigenvalues. 
(ICT) and (irm) . |f^.(o,w)r_ are also linearly independent for each given w 


Moreover, by 
G This 
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verifies the Kawashima condition (^) for the system (14.291) . Noting (14.211) . one can conclude 
the lemma from Theorem 1.1 of Shizuta-Kawashima [23|. □ 


5. Energy estimates with full dissipation 

This section is devoted to derive the energy estimates for the system (j3.8p which includes 
the full dissipative estimates of the component under the assumption (2IlDi) or ((Di). 
Throughout Sections [5H8] it will be assumed a priori that for sufficiently small 6 > 0, 

\\ u { t)\\He (5-1) 

5.1. Entropy dissipation estimates. One may first derive the energy estimates for the 
solution u by using the entropy fj = fi{u), defined by (|3.3U|) . 


Lemma 5.1. The following estimates hold, 
(i) Under i^i), 

d 
dt 

(m) Under (Di), 


J p{u) + < <5^ 


^ j ii{u) + \\u°\\l. 


~b 

2 


U 

L2 

U 


2 

L°° 


Proof. Applying {Vufi{u)y to (13.8p . the entropy assumption (Slz) implies 

d 


dtv{u)+ =Vufi{u) ■ Q{u). 


k=l 


Integrating (|5.4I) over by the entropy dissipation assumption ( 213 ), one has 


^ / p{u) + Ce Q{u) 


1,2 


< 0 . 


(5.2) 

(5.3) 

(5.4) 

(5.5) 


Since trivially 


2 2 

Q^{u) < Q{u) and recalling 0 dehned in (14.61) . it follows from (|5.5I) that 


dt 


7)(i2) + ||0°'u°||^2 ^ Q°(h)-0°M° 


L2 


(5.6) 


Now one may prove (i) under (QUSi). By the expansion ()4.1ip and the a priori assumption 
dSI]), using Sobolev’s inequality, it holds that 


Q^{u) - e^u° 


L2 


Y {u)vfu 


< 


L2 


IL2 


u 


<5^ 


u 


(5.7) 


Plugging the estimates (|5.7p into (|5.6p . by the assumption (2ti) (or equivalently, (I4.12p i. one 
deduces (|5.2p . 

0ne then turns to prove {ii) under ((Di). By the expansion (14.191) . it holds that 


Q^{u) — 0°h° 


< 


L2 


Y{u)vPu 


L2 


+ 


Y {u)ii!’v!’ 


L2 


< ll'W°lli2 + 


~b 

2 

~b 

U 

L2 

U 


Plugging the estimates ()5.8I) into (j5.6p . since 5 is small, one deduces (15.3p . 


(5.8) 

□ 


5.2. Energy estimates for One may next derive the higher-order energy estimates for the 
component vt by using the structure of the subsystem (13.261) . 

Lemma 5.2. For 1 < s < i, the following estimates hold. 

(i) Under i^i), 


^ [{A^uYA^’''{u)A^u^ + ||A*u°| 
dt J 


2 

L2 


< (5 




L2 


+ 


U 


+ 




(5.9) 
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(ii) Under (Di), 
dt 




2 

L2 


< (5 


( 


2 ,, ..,,9 

+ roo + 


') + 



v!’ 

V 


L2 II IH°“ 


L°° J 


1,2 



(5.10) 


Proof. For 1 < s < ^, applying to (I.S.26h and employing the commutator notation (IB.7I1 
yields 

d ^ d _ _ 

5i(A^h^) + = - J] [A^P>(n)] ui^+A^Q^u). (5.11) 

A:=l k=l 

Applying (A^'u^)'^A‘^’'’('u) to (15.lip , since by Lemmaand A^’^{u)A^’^{u) are symmet¬ 
ric, one has 

d 

dt (^{A^uYA^’\u)A^u''^+Y,{ i^"uYA^’\u)A’^’'’{u)A%^^ 


k=l 




= {A^iiY Idt (A0>(n)) + ^(^A^’\u)A'^’\u)'^ 


A^~> 


u 


k=l 


Xk 


(5.12) 


2Re (^(A"ix^)^A0’^(n) [A",A^'^(n)] n^J + 2Re ((A"n^)^A0’^(il)A"Q^(n)) . 


\k=l 


Here Ref = (/ + /^)/2. 

Now one may estimate the integration of the right hand side of (15.1211 over Replacing 
dtu by using the system (13.8p . sine £ > d/2 + 1 and by Sobolev’s inequality, the first term is 
easily bounded by 


L°o + l|Vu||^o 




u 


L2 


< (5 




2 

L2 


(5.13) 


By the commutator estimates (IB.81) of Lemma IB.31 using Cauchy’s inequality, one may bound 
the second term by 
d 


L2 


k=l 




< 

r\j 


A^ 


u 


L2 


x; (||v(/i*’^ti)) 


k=l 


L2 


L°° 




1,2 




A^{A’^^\u)) 


1,2 


Vn^ 


L°° 


(5.14) 


< (5 


A^it 


L2 


-h 


Vi2^ 


L°° 


For the last term, one can rewrite 

{A^u^fA^'\u)A^Q\u) 

= {A^itfA^'\u)Q^A^u^ + (A"u^)'^I°’^(il)A* {q\u) - &it 
= (A"i2^)'^I°’^(O)0^A"n^ + {A^i/^f [a^'\u) - I°’‘'(0)) 0^A"i2‘ 
+ {A^itf^A^'\u)A^ {q\u) - 0^n' 


(5.15) 


By (I4.23I) - (I4.24I) from Lemma HAl it holds that 

J 2Re ((A"{i^)'^I°’^(O)0^A"il^) = - J {A%^fM^A^u° < -Cm ||A^ 

0n the other hand, it is easy to bound 

/ 2Re ((A*i2^)^ (^A^’\u) - I°’^(0)) 0^A"n^) < ||it||^ 


-D||2 
Il2 


u 




L2 


< (i 




u 


2 

L2 


(5.16) 

(5.17) 
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Hence, integrating (I5.12p over it follows from (I5.13I) - (I5.17I1 that 


dt 


{A^uYA^’\u)A%^ + \\A%m, 


< 6 




u 


L2 


+ 




L°° 


+ 




u 


1,2 


A* {q\u) - 


(5.18) 


L2 


Now one may prove (i) under (2IJ2)i). By the expansions (|4.10p - (|4.11l) and using the product 
estimates ()B.9p of Lemma IB.31 one may estimate 


A" (Q\u)-Q^it 


< 


Y{u)u 


L°° 




< 5 


A^u^ 


1,2 


+ 


L2 

u' 


1,2 


A' 

+ 


u]u u 


A* (y{u)u 


1,2 


1,2 


u 


L°° 


(5.19) 


Plugging the estimates (|5.19p into (j5.18p yields (15.91) . 

One then turns to prove (ii) under (S)i). By the expansions (|4.18l) - (|4.19l) and using the 
product estimates ()B.9h of Lemma IB.31 one may estimate 


A" (q'’{u) - 
Y'{u)u 


< 


L2 
A^r^D 


U -\\^2 + 


A" [Y'{u)u^u 

A* fy'(i2)i2 


1,2 


+ 


o 11 r o 

L2 " 


A" (y'(u)u^u^ 

D 


L2 


+ 

< 6 


Y’{u)u'^ 


L°° 


A^id 


L2 


+ 




L2 


+ w 


11 ,°° 


+ 


A* {Y'{u)u^ 
2 

L°° 


1,2 


U 


(5.20) 


L°° 


U 


Plugging the estimates (|5.20p into (j5.18p yields (I5.10D . 


□ 


5.3. Recovering the full dissipation. Noting that Lemmas 15.1115.21 only contain the dissipa¬ 
tion estimates for the component vP, one may then use the partial Kawashima condition of the 
subsystem (I3.26p to recover the dissipation estimates for the component vP (and hence iP). 


Lemma 5.3. For 0 < s < i — 1, the following estimates hold, 
(i) Under (SUDi) and that d >3, 

_d 
dt 


—i 


w 


. T 

bl K 


d^ + 


^S+I~b 


U 


< \\Yu^l, + \\Y+^if\\l,+5 


(ii) Under (Di), 

A ^ 

dt 


—i 


w 


. T 

bl K 


A) d,e + 


< ||a^u°||2, + ||a^+V||^2+<5 




+ 


2 

L2 

ViP 


2 

L2 

ViP 


+ 


+ 


A^u^ 

2 



L2 



Here K = is defined in Lemma 

Proof. The subsystem (13.261) can be rewritten in the following form 
d d 

dtvf + j;i^Ao)4 = - E - A^Ao)) 4 + Q\u) := h\ 


(5.21) 


(5.22) 


k=l 


k=l 


(5.23) 
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Taking the Fourier transform yields 


+ i = i}. 


k=l 


Let K = K(^/\^\) be the skew-symmetric real matrix from Lemma 14.41 then 


(5.24) 




k=l 


k=l 


\k,b 


K >2\^\ {Ckin-I - diag(0, In-r)) ■ (5.25) 


Applying —i{u^)'^K to (I5.24D yields 

— Kdtv}' + {u^)^K ^kA^’'^{0)u^ = —i{u^)^Kh^. 

It then follows from (I5.25P that 


k=l 


- idt ({u^)'^KuA+ 2ck\^\ <2 1^1 




+ 2 




(5.26) 


(5.27) 


For 0 < s < £ — 1, multiplying (I5.27P by 1^1^^"*"^ and then integrating over the frequency space 
yields, by the Plancherel theorem, 


d 

df 




A' 


.+i^b 


L2 


< ||a*+V||^3 + 


A*+^i 


L2 




L2 


(5.28) 


To estimate ||A*/i^||^ 2 , one may use the product estimates (IB.91) of Lemma [6.31 to get 

A^((l^>(u)-I^>(0)) 4) 


< 


AX{u) - 


/c.b / 


L2 

A"Vu^ 


L2 


+ 


A" [a^'\u)-A^'^ 


L2 




(5.29) 


< 5 


A^+i 


u 


L2 


+ 




L°° 


On the other hand 


A' 


(<3‘(«)) 


L2 


< ||0°A*^'-° 


U 11^2 + 


A* (g^(i2) - 0^12^ 


L2 


(5.30) 


Hence, by (I5.29p - (l5.30p . it follows from (I5.28P that 

T 


df 




+ 


< || a ^ u °||;^2 + || a ^+ V ||22 + 5 






+ 


2 

L2 

A' 


.+i^b 


L2 


A* (Q^(i2) - 0^u^ 


(5.31) 


L2 


Now one may prove (i) under (22LDi) and that d > 3. By the expansions (I4.inp - (l4.1ip and 
the product estimates pB.9p of Lemma IB.31 it holds that 


A' {q\u) - 0^12^ 


< 


1,2 


A" 


I ft*’22 


< 


Y (22)ft 


< d 


A*+ift^ 


A"i2^ 


L2 


+ 


2d Y 


U 


(y{u) 

A" fy(ft)ft 


L2 


L2 


U 


(5.32) 


Plugging the estimates (15.321) into (I5.3ip . one deduces (I5.2ip . 
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One then turns to prove (ii) under (2)i). By the expansions (|4.18p - (|4.19l) and the product 
estimates (1B.9P of Lemma IB.31 it holds that 

A" (y'{u)u^u) + A" (y'{u)u^u^^ 




< 

rv-/ 


(q\u) - 
Y'{u)u 


< 

L2 ~ 

A^'U°||j^2 + 


L2 


L2 


AUY'{u)u 


AM(y'(u)-y' 


+ 


<5(||a^u°||^2 + ||u°| 


L2 


L°° 


)+5 


L°° 




+ 




+ 


Y\u)u^ 




L2 


AMy'(0)i2^ 


L2 


(5.33) 








L2 



Loo 


Plugging the estimates ()5.33p into (I5.3ip . one deduces ()5.22h . □ 

5.4. Synthesis. Now one can combine Lemmas l5.1H5.3l to arrive the following propositions. 

Proposition 5.4. For 1 < s < £ — 1, there exists an energy functional Ssiu^) equivalent to 
||A^u^||^^_^ such that the following estimates hold. 


(i) Under (2ir2)i) and that d > 3 

d 


+ + ||A^+iu" 


|2 


||A^n^||22 + 


u 


+ 




(5.34) 


{ii) Under (Di), 

<6(\\A^u^\\l, + \\u°\\l^ + 



2 

+ 


+ 



ht 


LOO 


L°° / 


L2 



(5.35) 


Proof. One may first prove {i) under (SUSi) and that d > 3 by using the (i) assertions of 
Lemmas l5.2H5.3l Taking first 1 < s < — 1 and then s = £vii the estimates ()5.9p of Lemma 15.21 
noting that ||/||j:/fe is equivalent to the sum ||/||j ;^2 + ||A^/||^ 2 ) one has 


^ J ((A*u‘')'^A0’^(i2)A"il^ + (AV)^A°’^(A)AV) + ||A"n°||^,_, 


< d 


A"y 


jje-s 


+ 


u 


+ 


Vtd 


(5.36) 


Next, taking first 1 < s < .£ — 1 and then s = — 1 in the estimates (15.211) of Lemma 15.31 one 

obtains 


_d 

dt 


J -i ((i^r ^'’Yk (icr+i + (ici'-i h^Y ^ (i^i' 


+ 


^s+l~b 


U 


2 


<\\A^uYHe-s+S 


u 


+ 


Vu' 


(5.37) 


For sufficiently small constant e > 0, if one defines the instant energy functional 
£^{u^) = j (^{A^uYA°’\u)A^u'' + {A^uYA^’\u)A^u^^ 


+ e -i 


(5.38) 


for 1 < s < £ — 1, since AP''’[u) is positive definite near u = 0 from Lemma 14.21 then it follows 
that £s{y^) is equivalent to ||A®u'’||^^_g. Moreover, since both e and 6 are small, one can deduce 
(ICTI) from (ICTP-(I07P. 

Similarly, to prove {ii) under (T>i), one needs only to use instead the {ii) assertions of Lemmas 
15.2115.31 in the arguments above to deduce (I5.35p . □ 
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Proposition 5.5. The following estimates hold, 
(i) Under (2ir2)i) and that d >3, 

rt 


Mmh + 

< 
rs_/ 

(a) Under (Si), 

Au^{t) 




(5.39) 


2 -U 
L2 + 


Au 


Ht-1 


u{mi2 + 


^ Jo 




< 

r\-' 


l^oll L2 + 



2 U / 


2 


2 

u\t) 

Auq 

+ / ( 

lt{T) 


u\t) 

+ 

Jo V 


L 2 

L°° 


(5.40) 


dr. 


Proof. One may first prove (i) under (QUSi) and that d > 3. Taking s = 1 in the estimates 
(I5.34P of Proposition 15.41 one obtains 

(5.41) 


dt 


< 6 


Aid 


L2 


+ 


U 


+ 




Taking s = 0 in the estimates (I5.2ip of Lemma 15.31 yields 


dt 


(^1^1 d.e + 


^ P L2 


| Af 2°||22 + 5 


Avt 

it 


2 

L2 


Vlt 


(5.42) 


Since I > d/2 + 1 and d > 3, by the interpolation estimates of Lemma lB.21 it holds that 


Vf2' 


< 

rsj 


Ait 


Ht-i 


(5.43) 


By this, integrating in time the sum of the estimates (j5.41h ~ (l5.42p and the entropy estimates 
(j5.2p of Lemma [5Tl since 6 is small, one can deduce (I5.39p . 

Similarly, to prove (ii) under (Si), one needs only to use instead the estimates (I5.35p of 
Proposition [521 the estimates (j5.22h of Lemma 15.31 and the estimates (15.3p of Lemma 15.11 to 
deduce (|5.4np . □ 


6. Full energy estimates 


In view of the energy estimates (l5.39l) - (|5.4Up of Proposition lH^l in Section^ it appears natural 
to derive the energy estimates for the derivatives of the remaining component u\. However, the 
evolution equation (I3.25P of ui is apparently not suitable for this goal since there would be a loss 
of derivatives if one directly performs the energy estimate. To get around this obstacle, despite 
that the component vf has already been controlled, one should still return back to derive the 
estimates for the derivatives of the whole u by using the original system p3.8p . 


Proposition 6.1. The following estimates hold, 
(i) Under (21LDi)-(2BS)2) and that d> 3, 

rt 




Ht + 


U 


+ 




+ 


A^-^ht 


1,2 


+ 


A^it 


L2 


dr. 


( 6 . 1 ) 


{ii) Under (2)i)-((D2) and that d> 2, 


Jo 


it 

T 11 it 11 r C.O T 

Vu^ 

+ 

> 

1 

+ 

A^it 


£,00 II 111 


Loo 


L2 



L2 


dr. 

( 6 . 2 ) 
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Proof. Thanks to the entropy estimates in Lemma FSTl by the Sobolev interpolation, it suffices to 
derive the energy estimates for A^u. Note that this is equivalent to derive the energy estimates of 
In the spirit of one dimensional study PIISIIITI, one may do the wave decomposition 

n 

Uxk = '^Wi^kfiiu,ek), k = l,...,d, (6.3) 

i=l 


where 

Wi^k = li{u,ek)uxt, i = k = l,...,d. (6.4) 

Using the product estimate (IB.9P of Lemma [6.31 one has 


A£-1 


u. 






L2 


^ ||wij,fc||^oo A^ ^(fi(u,efc) - fj(0,efc)) , + X] Ah{u, 


2=1 


2=1 


L2 


Sfc) II L° 


<d\\u\\H^ A^ 


2=1 


L2 


k = 1,... ,d. 


(6.5) 


It then reduces to derive the energy estimates for A^ ^Wi^k- But by (13.211) . it holds that 


Wj k ^k')^xk ^ ^ j 2,..., Ti, k 1,..., d. (6'6) 

j'=2 

This together with the estimates (I5.39p - (|5.4np of Proposition 15.51 implies that it suffices to 
derive the energy estimates of A^~^wi ki for each k' = 1,... ,d. Note that for the notational 
convenience the index notation k is changed into k'. 

Now it is crucial to derive the equation of wi^k'- Taking dx^, to the system (j3.8p yields 


d 

dt^Xy + 'y ^ (a {u)Ux^ = ^uQiu)Uxy ■ (6-'l') 

k=l 

Substituting the decomposition (16. 3p into (j6.7p , by (|3.17p and the system (|3.8p , one can deduce 

n n d 

Y.^tWkk'h{u,ek') +Y,Y1 Xi{u,ek) {Wi,k)xy 
2=1 2=1 k = l 

n 

= Wi^k'^uQ{uAi{u, efc/) 

2=1 


n d 


EE Aiu, ek)wi/^kri>{u, Bk) - Q{u) 

2=1 \2^ = 1 k = l / 

n d / ^ \ 

-EE Wi^kK{u,Bk)^ufi[u,ek) ( '^Wi',k'ri'{u,ek>) | 

2=1 k=l \2^ = 1 / 

n d 

-EE Wi,k (VuAj(u, Gfc) ■fii{u,ek') \ Wi/^k'ri{u,ek). 


( 6 . 8 ) 


i=l k=l 


Multiplying (16.8p by li{u) (owing to (5S)) and using (j3.18l) . one has 

d 

dtwi^k' + ^ Ai(u, Gfc) {wi^k)^^, 

k=l 

n n 

= ^/'i(u)VuQ(tt)fj(i2, efc/)u}j,fc/ - '^h{u)Vuri{u,ek>)Q{u)wi^k' 
2=1 2=1 

n d 

+ E E li{u)'^uri{u, ekf)ri/(u, ek)Xi'{u, Bk)wi,k'‘^v,k 

2,2^ = 1 k = l 


(6.9) 
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n d 


EE li{u)Vufi{u, ek)ri>{u, efc/)Aj(tt, ek)wi^kWi'^k' 




,k' •— ^ ^ Im- 


i,i'=l k=l 
n d 

- ^ ^ (ViiAi(t6, Gfc) • ri{u, efc/ 

i=l fc=l m=l 

One next needs to reexamine each term in (j6.9l) . Indeed, by (®), it holds that 

dQ 
' dui 


h = hiu)T^{u)wi^k' + y^Ji{u)VuQ{u)rj{u,ek')wj^k', 

1=2 


h = - li{u)Vurj{u, ek>)Q{u)wj^k', 
i=2 

n n d 

^=' = EEE !i{u)Vufj{u, ek')ri'{u, ek)Xi'{u, ek)wj^k'Wi',k-, 

j=2 i'=l k=l 
n n d 

^‘ = -EEE li{u)Vurj{u, ek)ri'{u, ek')Xj{u, ek)wj,kWi',k' 

j=2 i'=l k=l 


and 


d d 

uAi 


'd^l,k'^j,k' • 


( 6 . 10 ) 

( 6 . 11 ) 

( 6 . 12 ) 

(6.13) 

(6.14) 


h = — ^ ^k)wi^kWi^k' -EE( V^iAi(it, efc) • fj{u,ek' 

k=i j=2 k=l 

Lastly, notice that the left hand side of dlJl) is not a transported one, and it causes the 
inconvenience in doing the energy estimates. One may artificially recover this with the following: 


n 


{wi,k')^^ - {wi,k)^^, = [h{u)ux^,)^ - [h{u)u, 








X] (■“’ ®fc) ( '^uh{u) - (Vu!i{u)J j fi/(tt, ek')wi^kWi',k'■ 


(6.15) 


Moreover, for i = i' = 1 in the snmmation, by (®), 

fj{u,ek) (^uh{u) - (Vuh{u)^ ^ fi(u,efc/) = ej (vuli{u) - (Vuli{u)^ ^ ei = 0. (6.16) 

Hence, 

y (7^\ 

rf{u,ek)(vJi{u)-(Vuh{u)^ j ri>{u, ek')wi^kWi',k' (6.17) 
i+i'>3 X / 

Collecting these facts (I6.10p - (l6.14p and (I6.17P together with (16.61) . one may thus deduce from 
(j6.9l) that 

d 

dtwi^k' + ^ Ai(u, Gfc) {wi^k'). 


'^k 


k=l 


T f~\ i'~\ ~ \ ~ ~ X”'' o r~\~ 

h[u) — {u)wi^k' - 2 ^ -^{u,ek)wi^kWi^k' + 2 ^^ j,k'[u)Uj^x^, 


dui 

n n d 


k=l 


(6.18) 


i=2 


+ EE E Tji'^kk* 

j=2 i'=l k,k*=l 


Xk'^i',Xk* • 


Now one may prove (i) under (21LDi)-(2IJ£>2) and that d > 3. Assume (2IlS)i)-(2IlS)2), then 
it holds that ^-(hiei) = 0 and ||^(ttiei,efc) = 0. Using this fact, one may fix hi and view 
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^^{u) and Sfc) as functions of ■u^' = {u 2 ', • • •, UnY ■ By the mean value theorem, it holds 

that _ _ 

5Q, (,) ^ ^ ^ ViYu)uj = ^ Vij{i 


d\i , 


dui 


[UjUi 


i=2 


J=2 


and 


Also, 


^ ^ B/j-fc(u)uj = 2^ Vy,-fc(u)uj. 


dui 


i=2 


f=2 


Ai(u, efc) = Ai(Miei, e^) + ^ W*Yu)uj = Ai(0, e*,) + ^ 

i=2 i=2 

Hence, under (2IJDi)-(2D2)2), it holds that 

’ dui 


U)Ui 


hiu)j^{u)wi,k' = EE' iiiu)Vij{u)ujWi^k', 

i=l j=2 


and 


d n 


Wj^k{u)ujWi,kWi,k'- 

k=l ^ ^ k=l j=2 

Substituting (I6.22p - (|6.23p into (|6.18l) . one can get 

d 

dtwi^k' + Z 


(6.19) 

( 6 . 20 ) 

( 6 . 21 ) 

( 6 . 22 ) 

(6.23) 


'^k 


k=l 


n n n n d 

j=2 j=2 j=2i'=lk,k*=l 

The energy estimates on (j6.24p yields 

d 


(6.24) 


:=AAi. 


1 d_ 
2dt 




L2 


+ 


u, efc)^A^ A^ 


k=l 


(6.25) 




= / A^ W^A^ ^wii,fc'-Z / \Ai(u,efc) (u;i,fc')^^A^ 

fc=i ■ 

Integrating by parts, by (16.211) . one has 


j Xi{u,ek)(^A^ = \ j Ai(u, e^) 


A^ 


a;fc 


{W*iu)iYY, 


< 

r\j 


+ 




A^ 

2 


(6.26) 


m 


IIh- 


By the commutator estimates (IB.8h , the product estimates (IB.9h from Lemma rB.3[ the expansion 
(I6.2ip and since d > 3, £ > d/2 + 1 > 2, it yields 


A^ \Ai(u,efc) (uipfc') 


xk 


L2 


< 

r\_/ 


< 

/•N-/ 


< 


VAi(u,efc) A^ 

L°° 

^2+ ^(Ai(u,efc)) 

V (^lT*(u)u^) 

L°° 

A^ 

+ 

L2 

A^-i (^iy*(it)u^) 


2d 

L3=2 


2d 

iH=2 








L°° 


+ 




£,oo 


+ 




L2 




(6.27) 
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Lastly, one may apply the product estimates (IB.9P from Lemma lB.SI to obtain, by the expression 
of in (|6.24p , 



< 


+ 


+ 


+ 



L2 ~ 


L°° 


Loo 


L2 



(6.28) 


In view of these estimates (j6.26h - (l6.28p . it follows from (16 .25 11 that 


dt 




k' 


L2 


+ 




+ 




L2 


+ 




L2 


\u\\jje 


(6.29) 


Finally, integrating the inequality p6.29p directly in time, by ()6.3p - (l6.6p together with the 
estimates (j5.39p . one has 


^ \\uo\\h‘ + ( 


L°° 


+ 






+ 




L2 


+ 




u 


1,2 


) dr. (6.30) 


By Cauchy’s inequality, this implies p6.1D . 

One then turns to prove (m) under ((Di)~((D 2 ) and that d > 2. Assume ( 2 )i)~(!l) 2 ), then it 
holds that ^g^(h) = 0 and e^) = 0. By this fact, it follows from (|6.18p that 


dtwi^k' + ^ Ai(?2, efc) 
k=l 

= h{u) ^ {u)J{u)Q{u)wi^k' + ^ Sj^k'{u)uj, 


dui 

n n d 




i=2 


A ^ ^ Tji',kk*{u)Uj^xi,Uii^xi^* ■—Af . 

j=2 i'=l k,k* = l 


The A^-i energy estimates on p6.3ip yields 


iiL 

2dt 


^ k=i 

f f A^"\ Ai(it,efc) 


(6.31) 


(6.32) 


k'- 


k=l ■ 


Integrating by parts, by the linear degeneracy (T> 2 ), one has 



(6.33) 
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By the commutator estimates (IB. 81) from Lemma lB.31 by again ( 2 ) 2 ), it holds that 




L 2 


VAi(u,efc) ^ (wgfc') 

LCO ^ ' 

n 




1,2 


+ 


A^ Ai(A,efc 


2 Pd 




xk 




i=2 


L°° 


A^ 


fc' 


L 2 


(6.34) 


+ 5^||A^-2 (Vs,Ai(ti,efc)Vt 

i=2 


2Pd ||Vri)l,A;'||^Pd 




< 




+ 


A 


•-In' 


L 2 


+ 


A' 




L 2 


U 


j 


where = d for d > 3; = 4 for d = 2 and ^ > 5/2; = 2/(3 — /) for d = 2 and 2 < I < 5j2. 

Lastly, one may apply the product estimates ()B.9D from Lemma IB.31 to obtain, by the ex¬ 
pressions (I4.18p - (|4.19l) of Q, 

2 .. 


> 

1 

to 

< 



L2 ~ 




+ 





LOO 


L2 



+ « roo + 


In view of these estimates (j6.33l) - (l6.35p . it follows from (16.321) that 


L 2 


(6.35) 


d 

dt 




k' 


L2 


< 


u 


+ « roo + 



+ 

> 

1 

+ 



LOO 


L2 



L 2 


\u\\Ht ■ 
(6.36) 

This is an improvement of (I6.29P which then, together with the estimates (l5.4Up . leads to (16.2p 
in this case. □ 


7. ESTIMATES UNDER (D 3 ) 

This section is devoted to derive an L‘^ estimates of the solution in preparation to prove 
Theorem [3] under (2)i)-(2)3) and the initial condition uo,! G A*? with 1 < g < 2. However, one 
may not be able to derive the estimates of ui by directly using the subsystem (j3.25l) . The 
essential difficulty lies in the presence of the second summation involving vt’ on the left hand 
side. One feasible way is going back to the original system (13.81) of u and then resorting to the 
wave decomposition. 

In the spirit of the one dimensional study [aES], one may consider, owing to (*B), 

vi = li{u)u. (7-1) 

But by (13.211) . it holds that 

n 

Vi = Ui + '^ lij(u)uj. (7.2) 

i=2 

The key point is to derive the L'^ estimates of vi, and the result is stated as the following. 

Proposition 7.1. Under ( 2 >i)“((D 3 ) and that d >2 and vifi G L'^ with 1 < 9 < 2, it holds that 



Proof. First of all, one needs to derive the equation of vi. Multiplying the system (13.8p by li{u), 
by (I3.19p . (|3.2ip and an use of ()3.8p to substitute dtu, one can get 

d 

dm + '^Miu,ek){h)xk 

k=l 
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li{u)Q{u) + dtu^Vuh{u)u + ^ Ai(t6, ek)ul^Vuh{u)u 


k=l 




du, ^ 

j=2 i=l 
d n n 

+ EEE( Al('U, G/j.) Aj/(?/,G/^) j Wj'^f^Tijf(u^ 

k=l j,j'=‘2 *=1 


(7.4) 


dlij{u) _ 

an,; 


Here in the second equality the wave composition (16.31) for Ux^, is used. 

By ( 2 ) 3 ), (14.6p . (14.81) and (14.121) . it holds that 

^~ip(0) = 0) p = r + l,...,n, 

which implies that 

r n 

h{u)Q{u) = '^Iiq{u)Qq{u) + ^ “ ^~ip(0)^ Qp(n). 

q=l p=r+l 


(7.5) 


(7.6) 


Substituting (I7.6p and the expansions (I4.18l) - (l4.19p (due to (Si)) into (|7.4I) . using the expression 
dSSI) for Wji^k and (17.2p for ui, one may deduce 
d 

dtvi + Ai(n, efc)(ui)3 


/Xk 


k=l 

n 


d n 


Y1 Rip{^)nup+ Y, Rjp{u)ujUp + 

p=r+l p=r+l j=2 k=lj,j'=2 


'^j')xk (7-7) 


(UiUjUjiUjii 


ij'=2 j,j',j"=2 

The key point lying in (17.71) is that unlike (|3.25p there is no term involving v!’ on the left hand 
side, which makes it possible to derive the L‘^ estimate for vi. 

Now multiplying the equation by vi and then integrating over one has 


d II ~ nq ^ 


dt 


E 

k=l 


Xi{u,ek){\vi\‘^] 


Xk 


+ j (|f)ir|n°| + |nir"^|7i^|(|«°| + |Vu^|) + |ui|''|uf+ 

By the linear degeneracy (Sz), it holds that 


(7.8) 


Mu,ek)i\vi\'^)xk = X]^(«>efc)(uj)3,Jui|« 






L°° 


(7.9) 


By Holder’s inequality, the second integral in the right hand side of (|7.8p is easily bounded by 


Q II ~D|| I II ~ iia—1 

Li Ip IIloo + IPl Wli 


+ pi 


Li 


U 


U 

q-1 

Li 


L^i 


U 


“ llL2q + 




L^i 


L‘2i 


U 


(7.10) 


Plugging the estimates (I7.9p - (l7.inp into (17.81) yields 




^ iniiL'j 
A 


r + 




L^1 


u 


+ 


u 


L‘^1 


r Wl^i + 


Vi2^ 


L‘^1 


2 

L^1 


(7.11) 
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Applying the Gronwall lemma to (17.111) . one can deduce (|7.3I) . □ 

Remark 7.2. In the statement of Proposition pO] , it is required that G A*? with 1 < q < 2. 
However, this can be guaranteed under the assumptions of Theorem O Indeed, by (j7.2p and 
(ESI), it holds that 

r n 

Vl=Ul + '^ hq{u)Uq + ^ (llpiu) 

q=2 p=r+l 

which implies, since 1 < p < q <2, 

ll^'Rolli, < 11^1,olli, + ll^tolL, + ||«o||i2, < 11^1,olli, 

8. Decay estimates 

In order to close the estimates stated in Proposition 16.11 and Proposition 17.li this section is 
devoted to derive the decay estimates for the dissipative component vf. 


- lip{^)) Up, (7.12) 

+ II^oILp + II^oIIh^ ■ (7-13) 


8.1. Bootstrap estimates. In this subsection, one may derive the bootstrap energy estimates 
based on the family of scaled energy estimates with minimum fractional derivative counts derived 
in Section [SJ which implies that the decay of the higher-order norms of u^{t) can be deduced 
from the decay of the lower-order norms. 


Lemma 8.1. There exists a positive constant A > 0 such that the following estimates hold, 
(i) Under (21hDi) and that d >3, for 1 < s < d/2, 

2 


A^u^’it) 


< 




( 8 . 1 ) 


(m) Under ((Di), for 1 < s < i — 1 and 1 < s < d/2 + 1, 


A^u\t) 


^ —At 






Hl-s 


+ +<5||A"-i#(r)||5^2 +<5||u^(r)||^^^ dr. 


( 8 . 2 ) 


Proof. One may Hrst prove (i) under (2IlDi) and that d > 3. Recall the estimates (I5.34p of 
Proposition 15.41 and then Hx 1 < s < d/2 (and hence s < ^ — 1). Since i > d/2 1, 


u 


L°° 


+ 




< 


Loo 


A^iT 


Ht-s 


then by the smallness of <5, it follows from (I5.34p of Proposition 15.41 that 


-8s{u^) + ||A*n°||^,_, + ||A*+iu'^||;,_._, < <5 ||A^u' 


C||2 


lL2 • 


Adding ||A®u '^||^2 to both sides of (18.4p yields 


dt 


£s{u ) -|- 


A^ 


u 


Ht-s 


< llA^ 


-iC||2 


u 


lL2 • 


(8.3) 


(8.4) 


(8.5) 


( 8 . 6 ) 


Since £s{u^) is equivalent to it then holds that 

^^£s{u^) + X£s{u^)<\\A%%, 

for some constant A > 0. Applying the Gronwall lemma to (|8.6I) . one concludes (j8.ip . 

One then turns to prove (ii) under (Si). Recall the estimates (15.351) of Proposition 15.41 and 
then fix s so that 1 < s < ^ — 1 and 1 < s < d/2 + 1. Then it holds that 


ViT 


L°° 


< 


A^ 


u 


and 


-iDI| 2 < ||AS-lf,D||2 
LOO ^ ipr a Wnt- 


u 


Ht-s 
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It thus follows from (j5.35l) that 
d 




u 




<6i \\A^um, + \\A^-^u^\\l,+ 



+ 

A*ii^ 


v!’ 


L°° / 


L2 



Adding ||A®ii ‘^||^2 to both sides of (18.71) . by Cauchy’s inequality, one can obtain 


lL 2 

dt 


£siu ) + 




u 


ne-i 


< ||A^h"||"2+<5||A^-'n°||'2+5 


which implies 


d 


+ < llA^n'^ll' +<^||A 


s —1 ~D ||2 


dr""“ / ' / - II- “ iIl2 

Applying the Gronwall lemma to (18.9p yields (18.21) . 


11 


(8.7) 


( 8 . 8 ) 

(8.9) 

□ 


8.2. Linear decay estimates. In order to derive the decay of the lower-order norms of u^{t) 
appearing in Lemma l8.11 in this subsection one may consider the linear decay estimates for the 
linearized system of (13.26^ : 

( 8 . 10 ) 

[u’’ |t=o= 

where the linear operator £} is defined by 






( 8 . 11 ) 


k=l 


The solution operator e of the system (I8.10p has the following time-decay property: 

Lemma 8.2. Let d > 2. For l<p<2, a>0 and denote p* = rAm{2,dp/{d — p)}- There 
exists a constant A > 0 such that 


A“n'^e- 


Ur 


L 2 


<{l+t) Hp 27 2 


£> 

e “ Uq 


~b 


L 2 


_ d f 1 _ l\ _ a 

<{l+t) Hp 2 ) 2 


( 

n‘^^0 

+ 

LP 


LP* J 

+ e~^^ 




+ 

LP 


LPJ 

+ e-^^ 



L 2 


L 2 


and 




A"n°e-*^4 


L 2 


L 2 


_df i _1 \ _ a_ 1 

< (1 _g 2 Vp 2j 2 2 


_d( i_1\_a_i 
2{p 2j 2 2 


( 

n'^^o 

+ 

LP 


LP* ) 

+ e-^* 




+ 

LP 


LPJ 

)+e-"* 



L 2 


L 2 


where 11'^ and 11'^ are the projection operators defined by 

nV = (a2, ■ • •, flr)^, II°a'’ = {Or+I, ..., On)^, for a*' = ( 02 ,..., On)'^ G 


( 8 . 12 ) 

(8.13) 

(8.14) 

(8.15) 

(8.16) 


Proof. These decay estimates follow by refining the initial condition in [2]. Indeed, since similar 
as in Lemma 14.41 the linear system (|8.1Up satisfies the Kawashima condition, from Section 4.2 
of [2], one can decompose the solution u^{t) of (I8.10p as 

u^it) = e~^^^UQ = A'(t)'UQ -|- IC{t)uQ, (8-17) 

and the two components satisfy the following bounds: 
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and 




< Ce 


—ct 




( 8 . 20 ) 


for some constants C, c > 0. Then the estimates (I8.12p - ()8.15p follow in the standard way. □ 

8.3. Nonlinear estimates. In order to use the linear decay estimates of Lemma 18.21 one may 
rewrite the subsystem (j3.26h in the following perturbed operator form 

dtu^ + =M\ 


( 8 . 21 ) 


where the linear operator is defined by (18.lip and the nonlinear term 

d 


( 8 . 22 ) 


k=l 


For the convenience of the later analysis, one may record the following estimates of Af'^. 


Lemma 8.3. The following estimates hold. 

(i) Under (21bi)i) and that d> 3, for 1 < s < £ — 1, 


L 2 

for I < p < 2, 


+ 


As-i^b 


L 2 


< 5 




L 2 


A^ 


u 


L 2 


+ 




u 


Af^ 

<5{ 

A V p)^ 

+ 


+ 

it^ 

)+ 

A2(^ p)if' 


LP \ 


L2 


L2 


L°° / 



L 2 


for 1 < p < 2d/{d + 2), 


and 


Af'’ 

dp r 

A V Pj 

+ 

A V p) li*’ 


V 


L2 



L 2 




L 2 


< 6 




+ 


(m) Under (2)i), for 0 < s < i — 1, 


AW^' 

<(5 ( 
L2 ~ V 



+ 


L°° ( 


L 2 


+ + 




L 2 


+ 


U 




+ 


\L°° 


Af^ 

<5( 

A V p) 

+ 


)+ 

A V p)y'> 


it 


LP \ 


L2 


L'^J 


L2 



for 1 < p <2, 


and 


for 1 < p < 2d/{d + 2), 


Af^ 


L 2 


< 5 




\u 


\L°° 


+ 


U 


Af^ 

dp, ^ ^ ( 

.d(l-i)+2_k 
A V pJ +> 

+ 

A — -)+l ~D 

A V Pj y'^ 


A2(l p) + 2^b 


\ 


L2 


L^J 



2 

L 2 


and 

Af'^ 

for 2d/ {d + 2) <p <2, 


L 2 


< 5 




+ W 


) + 

it 


it 

/ 


L2 



AAf'’ 

<6 ( 
LP \ 


+ 

L2 

A'^ 


+ 

L2 


J 


+ 


1,2 

v!’ 

L°° ’ 





(8.23) 

(8.24) 

(8.25) 

(8.26) 

(8.27) 

(8.28) 

(8.29) 

(8.30) 

(8.31) 

(8.32) 
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and 


AAf^ 


L 2 


< (5 




L 2 


+ 




+ M 


\L°° 


+ 


u 


{in) Moreover, under (£>i) and that HuiU^^q < 6 with I < p < q < 2, for 1 < p < 2, 




J\f^ 


LP V 


LP 


< 

rs_/ 

+ 

< <5 

+ 


Vn' 


L 2 


+ w 




L 2 


L 2 

Il°° 


+ 


L^J 


Loo 


+ ll'W |li;,2 + 


+ 11^^ IIl2 + 


Aif 


Aid 


L2 


L 2 


and for 2 d/ (d + 2 ) < p <2, 


AAA^ 


< (5 


LP 




+ 




L2 


+ 




L2 


+ 


L2 


+ 




Hi 




L2 


L2 

, ifp < q; 


), ifp< q; 
), ifp = q, 


(8.33) 


(8.34) 


(8.35) 


LP 


< d 


A 2 


+ 1 l> 


L2 


+ 




+ « 


+ 


+ 




L2 


L°° 


+ 




m 


Avf 


L2 


L 2 / 
z/p = g. 


L 2 . 

(8.36) 

(8.37) 


Proof. One may first prove (i) under (2BSi) and that d > 3. In this case, recalling from the 
definition (j8.22l) of Ah'" and the expansion (I4.10l) - (l4.1ip . one has 

Af^ = Y {u){uVu^ + + u^u^) + Z{u)ituu, 


and 

Af'^ = Y{u){uVu^ + u^u). 

By the product estimates (IB.911 of Lemma fB.31 it holds that for 0 < s < £ — 1, 


(8.38) 

(8.39) 


Y{u)u 


AW^ 

L2 ■ - 

+ A*(y(?2)u 

< d (||a^+^^ 

and that for 1 < s < ^ — 1, since d > 3, 


A^+i 


u 


L2 


+ 


A^it 


L2 


VY 


A^vd 


+ 


u 


(8.40) 


L2 


+ 


vy 


L°° 


+ 


U 


L°° 


A^-i^^b 

< 

Y{u)u 

( 

A^fx^ 

2d -\- 



L2 ~ 


Ld- V 





+ 

< 6 


A^-\Y{u)u) 


( 


+ 

it 

L2 V 


Loo 



A^+i 


u 


L 2 


+ 


A^u^ 


L2 




2d 

L^ 


L°° 

+ 


(8.41) 


u 


For 1 < p < 2, by Holder’s and Sobolev’s inequalities, one can get 


Af^ 


LP 


Y {ii)u 


L2 


Vii 


2p A I j U II 2p 


+ 


< 

rsj 


< 

r\j 


u 


Y{u) 

2 

L 2 p 


U 


+ i( 


Vu' 


2 

L2 p 

b 


Z{u)u 


bli 2p 

L 2 "l^ 


U 


2p “h I ] I j 2p “h 


U 


A^V-i)it 


+ <5 


L2 


AA-i)u- 

+ 

A -')+! ~b 

A V py u 

+ 



L2 


L2 



(8.42) 
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For 1 <p < 2d/{d + 2) and hence dp/{d — p) < 2, d{l — 1/p) + 2 < d/2 + 1, it holds that 


Af^ 

dp A 

Y {u)u 

A 


dp + 

it 




L 2 V 


Ld-p-dp/2 



_dp_ 

Ld-p-dp/2 


< 6 


.dri-i')+2_k . dfl-iWl _k 

A V u'’ + A V pj^ u'" 


L 2 


L 2 


(8.43) 


Thus, (I8.4np - (l8.43h yield (|8.23p - (l8.25p correspondingly. On the other hand, (I8.26p can be proved 
with the simple modification in (I8.43p . 

Next, one may turn to prove (ii) under (Si). In this case, by the expansions ()4.18l) - (j4.19p . 
the nonlinear terms can be expressed as 

= Y{u){u'Vu^ + iPu) + Z{u)u^u^u 


and 

= Y{u){uVu^ + u^u + u^u^). 

By the product estimates (IB.9p of Lemma fB.31 it holds that for 0 < s < £ — 1, 




L 2 


< 


+ 

+ 

< 6 


Y {u)u 
k^{Y{u)u) 
Y{u) 




L 2 


+ A" A 




11,2 


L 2 




+ w 



A.^it 


vA 

+ 



L2 


LOO 


A%Yiu)-YiO)) 


L 2 


U 


1,2 


+ || a ^£(°||^2 + 




L°° 


+ \\u 


+ 




it 

+ 6 

it 


L2 


Loo 



2 

Loo 


For 1 < p < 2, by Holder’s and Sobolev’s inequalities, one can get 


Af^ 


LP 


< 

< 6 
< 6 


Y (u)u 


L 2 


Vu^ 






) + 

Z{u)u 


it 

2p 

it 

/ 


1,2 





Vit 

2p ^ 

it 

2p 

it 

) 


L'^-P ' L'^-P 




L°°/ 


, (i(l--')+l _k 

A V pj 

+ 


+ 



it 

) 


L 2 


L 2 


1,2 


L°° ) 


For 1 < p < 2d/{d + 2), it holds that 


Af^ 


ip_ ^ 




Y (u)u 


L 2 




dp ~t~ ^ ^P ' A 

L d—p—dp/2 L d—p—dp/2 


u 


2 

2dp 


< 5 


( 

. d( l-i)+ 2 _k 
A V p)^ 

+ 

. dfl-i^+l -D 
A V Pj^ u 


A 2 (l p) + 2 yb 

V 


1,2 


L^J 



1,2 


For 2d/{d + 2) < p < 2, by the product estimates ()B.9p of Lemma [6.31 it holds that 

IVn^i 


AAf'’ 

< 

Y{u)u 

dp 

V^it 

+ 

Y{u)u 


LP ~ 


Ld-P 


L<^ 



+ 

V{Y{u)u) 



Vlt 

+ 

Vit 

2p 


+ 

L 2 


L 2 
I ll~D|| 

2p “T rU 2p 

L'^-P ' L'^-P 


2p 


vy(u) 


L 2 


2p 

L 2 =P 


u 


< 6 

CN-/ 

+ <5 



+ 

. dri-i^+i _k 
A V p}^ u 

+ 



L 2 


L 2 



1,2 




L 2 


u 


Loo 


(8.44) 

(8.45) 


(8.46) 


(8.47) 


(8.48) 


(8.49) 


Thus, (|8.46p - (j8.49p yields (I8.27p - (l8.28h . (18.30^ . (I8.32p . respectively. On the other hand, (|8.29p . 
(j8.3ip . (I8.33P can be proved with the simple modifications in (I8.47p - (I8.49D . respectively. 
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Finally, one proves [iii) under (Si) and ||ui||j ;^5 < 5 with 1 < p < q < 2. By (|7.2I) . the 
expansions (I8.44p - (|8.45p can be rewritten as 


= Y ('u)i}i(V£(^ + u°) + Y{u)it{V^ + vP) + Z{u)ititu 


and 


= Y{u)vi{Vu^ + u°) + Y{u)u\VtP + u^). 
One may first assume p < q. For 1 < p < 2, it holds that 




LP 


Y{u)vi 


+ 


Vi2^ 


qp + tt qp 
Lq-p " "Lq-p 


< 


2p 

, 1,1 1 - 


Y (u)u 

s(^ 




L 2 


+ \u‘ 


+ 


L 2 


L 2 

Au^ 


+ 


L 2 

L^J 

+ 5 


+ 


Z{u)u 



2p 



L2 





L 2 


+ U 


lL2 


For 2d/ (d + 2) < p < 2, it holds that 
< V(Y(u)i 

LP 

Y {u)u 
Y{u 
VY{u 
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+ 

+ 

< 6 
r\j 
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vit 

qp 


qp 


Lq 



L9-P 

II "Lq-p 


dp 
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Y{u)u 

L2 1 

L' 

i-P 





IVU^II 2p 




L 2 


L°° 


2p 


Vi 2 ^ 




2p + 


1,2 


+ 


VA 

VY{u 


L 2 


U 


L 2 


2p 


2p 


U 


L°° 


(8.50) 

(8.51) 


(8.52) 


(8.53) 





+ 


+ 



L2 


L2 


L2 



L 2 




L 2 




L 2 


+ 


A ^ p'u 


L 2 




L 2 


+ 


u 


Here the following bound is used 

llVuilli, <||ui||^, + ||Vui||^. <d, (8.54) 

which follows from an use of Lemma [KTl Thus, (I8.52h - (l8.53p yield (I8.34h . (I8.36p . respectively. 
On the other hand, (I8.35h . (18.371) can be proved with the simple modifications in (I8.52l) - (|8.53l) . 
respectively. □ 


8.4. Nonlinear decay estimates. Now one can prove the decay estimates of ti^: 

Proposition 8.4. The following decay estimates hold. 

[i) Under (QUSi) and that d > 3 and 1 < p < 2, p < d/2, then for all s G [0,s|] O [0,d/2) 
with = d(l — 1 /p) + 1, 

||A^u'^(t)||^,_, < iLo(l + (8.55) 

and for all s G [ 0 , — 1 ] O [ 0 , d/2 — 1 ), 

\\A^u^it)\\^,_^<Koil + ty-^ip-^y^^-Y (8.56) 

(ii) Under (Si) and that d > 3 and 1 < p < 2, p < d; requiring further that p < 2d/(2(3 — 
1) + d) if i < 3, then for all s G [0, O [0, d/2 + 1) with = min{d(l — 1/p) + 1,1 — 1}, 

||A^u‘'(t)||j^,_. < Ko{l + (8.57) 

and for all s G [0, S 2 — 1] O [0, d/2), 

\\A^uyt)\\j^,_^ < Ko{l + t)--^^r-^y^^~Y (8.58) 
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{Hi) Under (£>i), (Ss) and that d > 2 and < S with 1 < p < q < 2, q < d; 

requiring further that p < 2d/{2{3 — i) + d) if i < 3, then for all s € [ 0 ,^ 3 ] n [0,d/2 + 1) with 
S 3 = min{d(l /2 + 1/q - l/p) + 1 , d{l - 1/p) + 2,£ - 1 }, 


\\A^u^{t)\\^,_,<Ko{l + t)--^^-r-"^^-K 

and for all s G [0, S 3 — 1] PI [0, d/2), 

||A*u°(t)||^,_, < Ko{l + 

Here Kq = ||uo||j|^f + ||'Wo||j;^p + ||^o|Lp* sufficiently small. 


(8.59) 


(8.60) 


Proof. One may first prove {i) under (SUSi) and that d > 3 and 1 < p < 2, p < d/2. Denote 
ap = d/2{l/p — 1 / 2 ) and take si = minls^, (i /2 — e} = min{(i(l — 1/p) + 1 , d/2 — e} < d/2 for a 
sufficiently small positive constant e. Then dehne 


Mi{t)= sup <{ (1 + 

0<T<t 


A^^u\t) 
2 


Ht-Si 




IL2 


+(1 + t)^-^ u^t) ^3 + (1 + t)^ 


’(^)l 


(8.61) 


First, taking s = si in (| 8 .ip yields 

0 

On the other hand, by the Duhamel principle, it follows from (I8.2ip that 

u\t) = + f e-(*-^)^A/’^(T) dr. 


(8.62) 


(8.63) 


Then one may apply the linear decay estimates (|8.12l) and (I8.14p of Lemma 18.21 to the formula 
(I8.63P to obtain 


1,2 


<(l+i) 2 , |P°|| * + 




L 2 


+ 




+ 






M^{r] 


LP 


+ 


M^{r) 


LP* 


+ 


(8.64) 


A"W(t) 


L 2 


+ 


A^l-l_^b(T) 


L 2 


dr. 


By the estimates (|8.23l) - (|8.26p from Lemma 18.31 it holds that 


<S 

< 6 


LP 


+ 




LP 
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L 2 


+ 


A^i- 1 ^^ 


L 2 




HI-H 




1 1,2 


+ 




Loo 


+ 


Loo 


+ 


A^p pJu^ 


2 

L 2 


(8.65) 


A"i 
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jjt-H II 11^^ 


+ 




2 

L 2 


Here in the first inequality it has been used that si < (i(1 — 1/p) +1 and that d{l — 1/p) + 2 < 
d/2 + 1 < i for 1 < p < 2d/{d + 2 ), and in the second inequality it has been used that, since 
Si < d/2 and I > d/2 + 1 and by the Sobolev interpolation, 


Vtd 


L°° 
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Loo 
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r\j 




*U 


Iji-n 


( 8 . 66 ) 
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Then plugging the estimates (|8.65l) into (I8.64p yields 


< (1 + + 






A^i 


u 


f^e-si II IIA2 


dr 


(8.67) 


+ 




+ t-T) “ 


A 2 


f i-ihA 


u (r) 


dr. 


L 2 


Now one may estimate the two time integrals in the right hand side of ()8.67h . In view of the 
definition of A4i(t), the first time integral is bounded by 

rt ~ ~ 

5 {1 + t- (1 + dr < 5(1 + ^/M^. ( 8 . 68 ) 

Jo 

Here it has been used that cr^ + si/2 > 1: if si = (i(l — 1/p) + 1, then ap + si/2 = d/4: + l/2 > 1 
since d > 3; if si = d/2 — e, then Up + si/2 = d/{2p) — e/2 > 1 since 1 < p < d/2 and one can 
take e sufficiently small. On the other hand, since d/2(l — 1/p) < si, by the interpolation the 
second time integral is bounded by 

rt 


Jo 




-2(7p-f (1- 


A4i(t) dr 


= /(l + t-r) ^{1 + t) 4A4i(r)(ir < (1+ t) 

Jo 

since si < d/2. Thus substituting the estimates (|8.68p - (l8.69p into (I8.67p yields 

||A^^'u'"(t)||^2 + 11^^^ ^ (1 + ^) ^ ^7^0 + 5y^ JAi{t) + A4i(f)^ . 

Then by (I8.7np . it follows from ()8.62p that 

A^^u\t) + j {K^ + + 7W?(r)) dr 

< (1 + t)-2-P-h (^2 ^ ^ _ 

Lastly and similarly, by using in addition the estimates (I8.26p . one can get 

||^'^(^)|lr2 ^ (1 + i) (II^oILp + II^oILp* + 


(8.69) 


(8.70) 


(8.71) 


L 2 

AA°(r) 


< 


T + (||AA"(r) 

rt ~ 

Jo 

(1 + (^Ko + 5^/M^ + Adi 




1,2 


dr 


(8.72) 


i(r) + Adi(r) dr 


and 


pt 

’(OIL 2 ^ (1+ + J (l + t-r)"'^J’"^(l + r)"‘^J’"^ (^5v^Adi(r) + Adi(r)^ dr 


< 


(1 + t) 2 ^iLo + 5v^Adi(t) + Adi 


(8.73) 


since si > 1. Now combining the estimates (I8.7UI) - (I8.73I) and recalling the definition of Adi(t) 
yields 

Miit)<Ki + 6^Mi{t)+Ml{t), (8.74) 

which implies, since both 6 and Adi(O) < Kq are small, 

Mi{t)<Kl (8.75) 


By again the definition of Adi(t) and by the interpolation, one can deduce (I8.55p - (j8.56h from 
(j8.75p since e > 0 can be taken arbitrarily small. 
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One then turns to prove (ii) under (3i) and that d > 3, l<p<2, p<d and further that 
p < 2d/{2{3 — £) + d) i < 3. Take §2 = min{s 2 , d/2 + 1 — e} for a sufficiently small positive 
constant e. Then define 


(8.76) 


Miit) = 1(1 + + ||A«-ifi”(T)||^d 

+ (1 + Tf-- ||ii‘(T)||(^ + (1 + + l ||«”(t)|P^,| . 

First, taking s = §2 in (|8.2p yields 

(l|A«ii'(T)||y + ||A*>-‘i“(T)||^. + |iA‘(U||(„) dr. 

(8.77) 

On the other hand, if 1 < p < 2d/{d + 2) and hence p* = dp/{d — p) <2, one may apply the 
linear decay estimates (I8.12p and (I8.14p of Lemma 18.21 to the formula (j8.63p to deduce 


(1 + i) 2 ( ||wq|Lp + ||uo|| _^ + 






L 2 


+ 




L 2 


(8.78) 


+ 



Af^ir) 

+ 

AA°(r) 

dp + 

A'W^r) 

+ 

AS 2 -i_^b(^) 

Jo V 


LP 


L'^ 


L 2 



L 2 


dr 


if 2d/{d + 2) < p < 2, one may apply instead the linear decay estimates (|8.13p and (I8.15P to 
deduce 


< (1 +1) “"p 2 + IIAuqII + A^^iio + A'^^ 


L 2 


+ 




+ t- 




LP 


+ 


AM^{t) 


LP 


+ 


(8.79) 


L 2 


A*W^(r) + A"2-W^(r) 


L 2 


dr. 


By the estimates (18.301) from Lemma [ 8 ^ it holds that for 1 < p < 2d/{d + 2), 

2 


J\f^ 


dp ^ S 


A^^u 


H^-^2 


+ 




ljt-S2/2 ’ 


(8.80) 


where it has been used that §2 < (1 — 1 /p) + 1 and that d(l — 1 /p) + 2 < d/2 + 1 < i ioi 

1 < p < 2d /(d + 2), and by the estimates (I8.32p - (|8.33l) that for 2d /(d + 2) < p < 2, 


AN^ 


LP 


< 5 


A^^it 


Ht-S2 II I To 


) + 




/ 


L 2 



(8.81) 


These two estimates together with the estimates ()8.27p - (l8.29p and ()8.33l) imply that the two 
sums inside the two time integrals in the right hand side of ()8.78l) and ()8.79p are bounded by 


A^'b 


+ 

< 5 


b 

'2 -b 


^ _ +||A^'^-V||^,+ 

Ui-S2 II iir^ 

2 




+ W 


A 2 li 


j^t-S2/2 


+ 


u 


A^ 2 - 




H^-^2 




III 



+ 

L2 


L°°) 




)- 

A ^ 'b 

A 2 u 

2 

-* 2/2 

it 

L°° 

A52-l^b 


(8.82) 


jje-S2+i 


Here in the first bound it has been used that §2 < d(l — l/p) + l, §2 < d/2 + 1, and in the 
second bound it has been used that, since §2 < d/2 + 1 < i and by Sobolev interpolation. 




L°° 


+ \\u 


< 

Il°° 


A^^it 


h^-h 


U S2-1 ~D I 


U 


L 2 • 


(8.83) 
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Plugging the estimates (|8.82p into (18.781) or (I8.79p . in each case one can get 




L 2 


< (1 + + [\l + t - 

Jo 


dr 


A ^2 


u 


fjt-S2 II IIA2 


dr 


+ 

pt 

s 

A- 

2 ~b 
u 

J 

0 

rt 




0 


LOO 


(8.84) 


AS2-i£^b 


fjl-S2 + l 


dr. 


Now one may estimate the three time integrals in the right hand side of (I8.84p . In view of the 
definition of A 42 (t), the hrst time integral is bounded by 


<5 [\l + t - (1 + dr < <5(1 + /M^. (8.85) 

Jo 


Here it has been used that cjp + S 2/2 > 1: if S 2 = c?(l — 1/p) + 1, then ap + S2/2 = (i/4 + 1/2 > 1 
since (i > 3; if S 2 = d/2 + 1 — e, then Up + S 2/2 = d/{2p) + 1/2 — e /2 > 1 since 1 <p < d and e 
can be taken sufficiently small; if S 2 = ^ — Ij then Up + S 2/2 = d/2(\/p — 1 / 2 ) + (£ — l )/2 > 1 , 
which trivially holds when ^ > 3, and follows from 1 < p < 2(i/(2(3 — i) + d) when i < 3. On 
the other hand, by the interpolation the second time integral is bounded by 

(1 + i - + r)-2-^-^A42(r) dr < (1 + M 2 {t). (8.86) 

To bound the third time integral one may use the interpolation to estimate 




< 

Amin{|-e,S2} ~b 


LOO ^ 



H 


! —min{ A—e,S 2 } 


< 


(1 + 1 ) 


min{^-e,g 2 } 




(8.87) 


Hence, the third time integral is bounded by 



+ i - r) 


.£2 


(1 + r) 


-2(7T) - 


min{^—e,S2} + S2“l 


M.2{T)dT 


<(l+t)-p-f>l 2 (i). 


( 8 . 88 ) 


Here it has been used that §2 > 1 and 2cjp + ((i/2 + S 2 —1)/2 = ap + S2/2 + {d/p — l)/2 > ap + S2/2 
since 1 < p < d and one can take e sufficiently small. Thus substituting the estimates (|8.85p . 
(j 8 . 86 l) and (| 8 . 88 l) into (18.841) yields 


|A^^'u*"(t)||^2 + ||A^^ ^'^*^( 011 x 2^(1 + ^) ^ (^Ko + J^2{t) + A42(t)^ 


(8.89) 


Then by (18.891) together with an use of the estimates (| 8 . 88 p , the bound of the third time integral 
in the right hand side of p8.84p . one may deduce from (I8.77P that 

(kS + + Ali(*)) dr 


(1 + t)-2-P-*2 ^ A ^ 


Finally, note that one needs to revise the estimates p8.72p - (l8.73p in the case (i) since now 
(Tp — Si > 1 may not be true for <7/2 < p < <7. Then it is crucial to exploit the exponential 
decay property of the high frequency part in the linear decay estimates of Lemma 18.21 Indeed, 
applying the linear decay estimates p8.12l) and p8.15p and using in addition the estimates p8.3ip , 
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one can deduce 

~ (1 + ^) (II^oIIlp + II'^oILp* + 


u 


L 2 


Ar°(r)||^^^) ||AA^(r) 


< {1 + Kq + J { 1 + t - Ty^{l+ Ty ^~'2 (^ 6 ^/MyT) + M 2 {t)^ dr (8.91) 

+ re-^(‘-")(l + r)- 

Jo 

< (1 + (^Ko + 6 s/M 2 {t) + M 2 (t)^ , forl<p< 


min{a-e,S 2 } 

- 2 - M2(T)dT 


2 d 

d + 2 ’ 


|^^(^)||2,2 ~ (1 + ^) (II^oILp + II^^oIIlp + "“0 ^ 2 ) 


LP 


+ 


+ ^ (l + t-r)-^^ (||AA‘'(r) 

< (1 + t) (^Kq + syM- 2 (t) + M -2 


AAA°(r)||^J ||AA^(r)||^2 (8-92) 


, 2 d 


and 


p't 

’(OIL 2 < ( 1 + + J (l + i-r)"'^J’"^(l + r)"‘^J’"^ + Ad 2 (r)^ dr 


+ / e-^(*-")(l + r) 


— 2(7T) — 


min{^-e,S 2 } 


2 Ad2(r)dr 


(8.93) 


< 


(1 + i) 2 ^ii'o + <5\/A^2(i) + -^2 


Here it has been used that S 2>1 and that ap + d/4 = d/{2p) > 1/2 since 1 < p < d. Based on 
these estimates (I8.89p - (j8.93p . one can get that for d and A 42 ( 0 ) < Kq small, 


M 2 {t) < Kl 


(8.94) 


This implies (I8.57p - (l8.58h . 

Finally, one turns to prove {iii) under (T>i), (TI 3 ) and that d > 2 and ||hi||^q < 6 with 1 < 
P < q <‘2, q < d and further that p < 2d/ (2(3 — £) + d) if < 3. Take S 3 = min{s 3 , d/2 + 1 — e} 
for a sufficiently small positive constant e. Then define 


M^{t)= sup l{i + rf-p+^^( \ , +||A^'3-V(t)||' 1 

o<T<t L V “3 ^ y 


Ml + rf” i\r) „ + (l + A 


’Wll 


L 2 


(8.95) 


One can still use the same estimates (I8.77p - (|8.79l) in the case (fz), however, the nonlinear 
estimates can be improved by using the estimates (|8.34I) - (I8.37I) from Lemma 18^ Indeed, by 
the estimates (|8.34p - (j8.35p . it holds that for 1 < p < 2, 




LP 


<5 




II IIL 2 


+ 




L 2 


U 


+ 


Ay 


L 2 


+ rd 


Il 2 


(8.96) 


and by the estimates (I8.36p - ()8.37p . it holds that for 2d/(d + 2) < p < 2, 


AAA° 


LP 


< 5 


A^H 


_ + ||A*'3-i^^d|| 

II Iin 


+ 


. d(l-i) ~t> 
A ' p'u 


L2 


L°° 


+ 




m 


Ay 


(8.97) 


L2 
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Hence, noting further (18.291) and (I8.33p . one may improve the estimates (I8.82p to be 




+ 


+ 




A^3 




A^*-^ p)+2 


2 


+ 


1,2 


+ 




m 


Aft^ 


L 2 


(8.98) 


1,2 


+ ft 


11,2 


Plugging the estimates (|8.98p to p8.78p or (I8.79p . in each case it holds that 

'-5 


< (1 + t)—p-^i;^o+ [ (1 + t - 

Jo 


+ 




10 


u\t) 




Ht-s^+l 


+ 




L 2 


(8.99) 


+ 


A 2 




u\t) A'^^^ p^i2^(r) ( Au\t) + \\u^{t)\\^ 2 )^ dr. 


Now one may estimate the two time integrals in the right hand side of (I8.99p . In view of the 
definition of the first time integral is bounded by 

+ t - r)-‘"^-^(l + r)-'^p-^VM 3 (r)dr < <5(1 + 




( 8 . 100 ) 


Here it has been used that ap + § 3/2 > 1: if S 3 = d{l/2 + 1/q — 1/p) + 1, then ap + S 3/2 = 
d/{2q) + 1/2 > 1 since 1 < 5 < <i; if S 3 = d(l — 1/p) + 2, then ap + S 3/2 = d/4 + 1 > 1; if 
S 3 = d/2 + 1 — e, then ap + s^/2 = d/ ( 2 p) + 1/2 — e /2 > 1 since 1 < p < d and one can take e 
sufficiently small; if S 3 = — 1, then ap + S 3/2 = d/2(l/p — 1/2) + (^ — l)/2 > 1 which trivially 

holds when i > 3, and follows from 1 < p < 2d/(2(3 — i) + d) when i < 3. To estimate the 
second time integral, it is noted that, similar to (I 8 . 88 p the first term is bounded by 

(l + t)-'"p-f A43(t). ( 8 . 101 ) 

Since d(l — 1/p) < S 3 one may use the interpolation, together with the estimates similar to 
P8.87I) . to bound the second term by 


/' 


(l + t-r)-'"p"^(l + r) 




dr 


( 8 . 102 ) 




Here it has been used that 2ap + d/2(l — 1/p) + d/4 = <Tp + d/2 > ap + {d/2 + l)/2 > ap + S 3/2 
since d > 2 and S 3 < d/2 + 1. Similarly, the last three terms are bounded by 


I 

Jo 


(1 + t —r) ""p 2 ^ (1 +1) p) ^A43(r)dr 


f 


= I (l + t-r) 2®(l + i) 4 2A43(T)dr < (1+ t) '’2 


(8.103) 


since S 3 < d/2 + 1. Thus substituting the estimates (|8.100p - (l8.103l) into (18.991) yields 

||A'^®?i'^(t )||^2 + ^^°(^)|| 2,2 ~ (1 + 0 -Aioit) + ■Aioi't)'^ ■ (8.104) 

Then by (18.1041) . similarly as (I8.90D . it follows from (I8.77D that 


A^^u\t) 




< (1 + {Kl + d^Mait) + Ml(t)) 


Finally, similarly as (I8.9ip - (l8.93p in the case (ii), one can get 

'^"(t)IL 2 < (1 +(i^o + <5/M^ + Ms 


u 


(8.105) 

(8.106) 
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and 

\Fit)\\L2 < (1 + {Ko + + Msit)) . 

Thus combining (I8.104p - (|8.107l) . one can get that for 6 and A^sPO) < Kq small, 

Ms{t) < Kl 

This implies (I8.59|) - (I8.6UI) . 

Remark 8.5. In the statement of the assertions (ii) and {in) of Proposition \8.4[ 
used implicitly that for 2 d/{d + 2 ) < p < 2 , 

II^^^oILp ~ II“o|L 2 + ||^■^o|L2 < A'o, 

which follows from an use of Lemma \B.1[ 


(8.107) 

(8.108) 
□ 

it has been 
(8.109) 


9. Proof of the theorems 


This section is devoted to prove the main theorems. By the standard local existence theory 
[mds] and a continuity argument one needs only to show how to close the a priori estimates. 

Proof of Theorem [ij Under the assumptions of Theorem[Tl one assumes a priori that ||'u(t) || < 

6 . Then by the (i) assertions of Propositions 16.11 and 18.41 taking si be the one in the proof of 
the assertion (z) of Proposition 18.41 it holds that 

rt 


\u{t)\\He ^ 


I'^olljyZ + 


A"iu'(r) 




dr 


(9.1) 


< 


\uo\\fjt+Ko / (l + r) 2 dr<ii:o. 


This closes the a priori estimates by taking Kq sufficiently small. Consequently, by the estimates 
(19.11) and the (i) assertions of Propositions 15.51 and 18.41 one can obtain (I2.8l) - (|2.10p and thus 
conclude the proof of Theorem [TJ □ 


Hi < 


Proof of Theorem [H Under the assumptions of Theorem[2j one assumes a priori that \\u 
5. Then by the {ii) assertions of Propositions 16.11 and 18.41 taking S 2 be the one in the proof of 
the assertion [ii) of Proposition 18.41 it holds that 

rt 

\h^ 




+ I + 11^'' '^°W|L2 + ||^^^)||'^) dr 


< 


+ ATq 


f 


(9.2) 


1 + t)"'"p"^ + (1 + r)" 2 ‘"^" 2 +") dr < a:. 


since d > 3. This closes the a priori estimates by taking Kq sufficiently small. Consequently, 
by the estimates (19.2p . the assertion (i) of Proposition 15.51 and the assertion {ii) of Proposition 
one can obtain (I2.13p - (|2.15p and thus conclude the proof of Theorem [2l □ 


Proof of Theorem 0 Under the assumptions of Theorem[3l one assumes a priori that ||h(t) || + 

||hi(t )||^9 < 6 . Then by the assertion {ii) of Proposition [ 6 T] and the assertion {in) of Proposition 
18.41 taking S 3 be the one in the proof of the assertion {Hi) of Proposition 18.41 it holds that 


u 




+ I + 11^'" '«°W|L2 + ||^^('^)||^^) dr 


< 


+ Kq 


f 


l + r)"'"p"^ + (l + r) 


(9.3) 




dr < K, 


0 


since 2ap + dj2 = d/p > 1 for 1 <p < d and one can take e sufficiently small. Moreover, by the 
estimates (j7.3p of Proposition 17.11 and Remark 17.21 it holds that 


II^i(^)IIl 9 ^ exp(CR:o) X Ro + 


f 










+ 


U 


<Ko + K^ 




2\p q) 2\p q 


dr < K< 


2 

L^q 

0 


dr 


(9.4) 
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since (i/2(2/p — 1/g) + 1/2 > d/ {2p) + 1/2 > 1 and d/2(3/p — 1/q) > d/p > 1 iov 1 < p < q < d. 
These close the a priori estimates by taking Kq sufficiently small. Consequently, by the estimates 
(insi), the assertion [ii) of Proposition 15.51 and the assertion [in) of Proposition 18.41 one can 
obtain (|2.18I) - (I2.20I) and thus conclude the proof of Theorem [3 □ 


10. Application to Damping Full Euler System 


In this section, we verify that the damping full Euler system of adiabatic flow satisfies all 
the assumptions (2li)-(2l4), (55) and (1Di)-(1D3). The system can be written in terms of the 
density /?, the velocity v and the entropy S of the flow as 

dtp + div(/??;) = 0, 

< dt{pv) + d\Y{pv +pld) =-pv, (10-1) 

dt{pS) + div(/?i;S') = 0. 


It is well known that the thermodynamics variables p, p, S, the temperature 6 and the internal 
energy e satisfy the relation 

de = edS-pd{-), (10.2) 

P 

and they can be determined by knowing any two of them; we may view p, 6, and e as functions 
of p and S. It is assumed that e is strictly convex with respect to p and S. 

The system (llO.ip can be rewritten into 

d 

dtS + '^VkSx^^ = 0 , 

k=l 

d d 

^ dtp + '^VkPxt,+'^p{vk)xk =0, (10.3) 

k=l k=l 

d 

dtVi + —Pxi + ^ Vk{vi)xf, + —Sxi = -Vi, i = 1,..., d. 

^ k=i ^ 

Denote u = (S', p, By a direct calculation, the matrix 


fv ■ (jj 0 \ 

A{u,oj) = 0 V -oj puj'^ I 

\t^ v-^Id] 

possesses d + 2 real eigenvalues 

\i{u,u) = ■ ■ ■ = \d{u,uj) = V ■ u, \d+i{u,i^) = V ■ UJ + y/p/, \d+ 2 {u,oj) = V ■ u - y/p/ 

and a complete set of left eigenvectors 

li{u,u) = (1,0,0^), lk{u,u) = (0,0,Xfc) (A: = 2,...,d), 

T) T) 

ld+l{u,Uj) = {—=,y/p/,pu}), ld+2{u,Uj) = i—=,y/p/,-pU}) 

y/Pp y/Pp 

and right eigenvectors 

I , rfc(u,a;) = j 0 I (/c = 2 ,... , d), 

0 / \XkJ 


ri{u,uj) = 


rd+i{u,uj) 



rd+2{u,u]) 



Here Xk (A; = 2,..., d) are d — 1 orthogonal unit vectors that are perpendicular to u. 

It is clear that the assumption (IB) is satisfied. Noting that Q{u) = (0,0,—u^)^, we check 
that around the equilibrium state rt* = (S*,/9*,0^)^ with S* and p* > 0 given constants. 
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• Take the entropy fj = + pe{p, S) and the entropy flux = v^{^p\v\‘^ + pe{p, S) + 

p{p, S)), and then dehne 

= 7:P\v\'^ + pe{p,S) - pe{p^,S^) - pS6{p^,S^) + p^S^6{p^, S^) - —p(/0*,5*) +p(p*,5*). 

2 

• Vuviu) • Q{u) = -plup < -^|up. 

• ^ j 3 — 

• ^uQ{u)ri{u) = 0 and V„Ai(u,w) • ri(u) = 0. 

• li{u)V uQ{u) = 0. 

These verify the assumptions ( 2 li)-( 2 l 4 ) and ( 1 i)i)-(Ti 3 ). Thus, Theorem [3] produces the unique 
global solution to the damping full Euler system in d> 2. 


Appendix A. Linear transformation 

In this appendix it will be showed that for the system a under the assumptions (2li)~ 
( 2 I 4 ) and (53), the conditions (|3.ip and (|3.2p can be always satisfied after a proper linear 
transformation. 

By (2li), one can apply to the system (II.ip the following linear transformation 


u = u{u*) 



Ir 

-1 

r+l<p,p'<n 



l<q<r 

r-\-l<p<n 



(A.l) 


It is easy to check that the system after this transformation keeps the same form as dni), and 
all the assumptions (2ti)-(2l4) and (53) still hold. Moreover, one has 


du* 


(0) =0, g = I,...,r; i = l,...,n. 


(A.2) 


This verihes the condition (13.11) . 

Now, by (2li), ( 214 ) and (13.ip . it is direct to get 


ri(0) G span{ei,... ,e,.}, 


(A.3) 


which implies that r*;^( 0 ) = 0 (p = r + 1 ,..., n) and for some A £ { 1 , • • •, 

<i(0)/0. (A.4) 

Then one may choose {z 2 , ■ ■ ■, = {Ij ■ ■ ■ t} \ {* 1 } and have that 

span{r)|'( 0 ),ei 2 ,... , 6 *^, 6 ^+!,... ,e„} = M”. (A.5) 

One can further apply to the system ([HD the linear transformation 

u* = u*{u**) := (ri( 0 ),ej 2 ,.. ., 6 *^, 6 ^+!,... , 6 ^) u**. (A. 6 ) 

Again, it is easy to check that the transformed system keeps the same form as and all the 

assumptions ( 2 ti)~( 2 l 4 ), (53) and the condition (13.11) still hold. Moreover, one has 


rr(0) =ei. 


(A.7) 


This verifies the condition (13.21) . 
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Appendix B. Analytic tools 

First recall the classical Sobolev interpolation of the Gagliardo-Nirenberg inequality. 
Lemma B.l. Let 1 < ( 7 , r < oo and 0 < a < /?. Then it holds that 


I|a“/IIl. < 


where 


for all 6 in the interval 


1 a 
p d 


A^/ 

/3 


i-e 

LI ’ 


1 

r d 


+ ( 1-4 


(B.l) 

(B.2) 


<e <i 


m 


with the following exceptional cases. 

(1) If a = 0, r/d < d, q = 00 , then one needs an additional assumption that either f tends 
to zero at infinity or f G L‘^ for some finite q > 0. 

(2) // 1 < r < 00 , meanwhile fi — a — d/r is a non negative integer, then (jB.ip holds only 
for 6 satisfying a/fi < 6 < 1. 

Proof. See THEOREM of Lecture II in [2l]. 


□ 

In the lemma above the restriction 0 < a < /3 is inconvenient for the use in this paper, one 
may employ the following variant. 


Lemma B.2. Let 2 < p < 00 and a, (5,^ S 

l|A“/||i. < 


Then it holds that 


e 


A^/ IIAVIll 


-0 

L 2 


Here 0 < 9 < 1 (if p = 00 , then it is required that 0 < 6 < 1) and a satisfies 

= (36 + 7(1 - 6). 


,, 1 1 
a + d I - - 


.2 P. 

Proof. Eor the case 2 < p < 00 , it follows from the classical Sobolev inequality that 

A^/ 


l|A“/||p< 


. , . ,, 1 1 

with Q = a + d { - — 


P 


By the Parseval theorem and Holder’s inequality, we have 


A7 < A»f 


IIAYIIP' 


where 0 < 0 < 1 is defined by (IB.4p . Hence, (|B.3p follows by (IB.5li - (IB.6li . 
Eor the case p = 00 , one may refer to Exercise 6.1.2 in [ 6 ]. 

Recall lastly the following commutator and product estimates: 

Lemma B.3. Let I > 0 and define the commutator 


(B.3) 

(B.4) 

(B.5) 

(B. 6 ) 

□ 


Then it holds that 


A\9 


In addition, for I > 0, 


A\gh) 


LPO 


LPO 


A\g h = 

A^(ci/i) -gA^h. 

(B.7) 

< \\^9\\lpi 

A^~^h + A^g ||h||^p4. 

LP2 LP3 ^ 

(B. 8 ) 

^ 5 LPI 

h}h + h}g h . 

LP2 LP3 ^ 

(B.9) 


Here po,p 2 ,P 3 G (l,oo) and 


1 _ 1 1 _ 1 1 

PO Pi P2 P3 Pi' 


Proof. One may refer to Lemma 3.1 in m- 


□ 
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